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THEORY  OP  NOISE  IN  A  MULTIDIMENSIONAL  SEMICONDUCTOR 
WITH  A  P-N  JUNCTION 


Prepared  by: 
Max  Solo* 


ABSTRACT:  This  thesis  discusses  the  fluctuations  of  noise 
in  a  two  and  three  dimensional  semiconductor  containing  a 
p-n  Junction.  We  consider  a  rectangular  parallelepiped 
single  crystal.  It  la  bisected  in  the  longest  dimension  by 
a  p-n  Junction.  Since  this  dimension  is  several  diffusion 
lengths  It  can  be  considered  infinite.  In  the  transverse 
plane  we  investigate  the  case  where  both  dimensions  are 
finite,  and  then  the  case  where  one  is  finite  and  the  other 
Infinite.  In  the  p-n  Junction  the  noise  is  the  result  of 
fluctuations  In  the  minority  carrier  density.  In  a  p-n 
Junction  there  are  two  classes  of  minority  carriers:  1. 
holes  in  the  n-type  material,  2.  electrons  in  the  p-type 
material.  Since  both  hole  and  electron  density  fluctuations 
are  similar,  we  discuss  only  the  former  In  detail.  We 
Investigate  the  differential  equations  for  a  two  and  three 
dimensional  semiconductor  with  a  p-n  Junction  and  find  the 
inhomogeneotis  form  of  these  equations.  These  equations  are 
solved  with  the  help  of  the  scalar  and  tensor  Green's 
function.  The  noise  problem  is  solved  by  using  these 
equations  as  Langevin  equations  and  interpreting  the  dis¬ 
tributed  sources  as  random  forces.  Then  the  noise  current 
spectrum  is  determined  with  stochastic  process  theory  after 
deriving  the  sources  from  basic  physical  models  and  the 
theory  of  stationary,  ergodlc,  Markovian  processes.  We 
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consider  two  cases  of  surface  recombination  velocity  on  the 
transverse  surfaces:  infinite  s  and  finite  s.  For  the 
infinite  case,  we  get  the  exact  solution  which  provides  an 
upper  bound  for  the  noise  spectrum  for  large  s.  For  an 
arbitrary  s  we  get  a  solution  but  have  confidence  in  the 
solution  for  only  small  s. 

Therefore  we  have  obtained  a  complete  solution  for  the  two 
cases  of  practical  interest:  large  and  small  surface 
recombination  velocity.  These  cases  should  prove  of  interest 
In  the  analysis  of  noise  phenomena  in  semiconductors. 
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This  report  describes  a  theoretical  study  of  the  fluctu¬ 
ation  noise  from  a  two  and  three  dimensional  semiconductor 
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problem  with  partial  support  from  PR-21,  The  report  13 
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CHAPTER  I 


INTRODUCTION 


Semiconductor  noise  studies  give  useful 
information  about  the  basic  physical  micros¬ 
copic  processes  in  semiconductors  and  in  the 
solid  state.  Furthermore,  noise  becomes  very 
inportant  when  a  semiconductor  device  such  as 
a  transistor  is  used  with  signal  levels  com¬ 
parable  to  the  noise. 

A  semiconductor  crystal  which  is  p-type 
at  one  end  and  n-t.ype  at  the  other  has  a 
transition  zone  which  is  called  a  p-n  junction 
(Shockley1,  Kittel^).  Current  is  carried 
across  the  Junction  by  minority  carriers;  that 
is,  electrons  in  the  p-type  region  and  holes 
in  the  n-type  region. 

Petrltz^jH  has  shown  that  noise  in  a 
p-n  Junction  arises  from  fluctuations  in  the 
concentration  of  minority  carriers.  Consid¬ 
ering  a  p-n  junction  as  an  ideal  one  dimen¬ 
sional  structure,  he  has  derived  expressions 
for  this  noise.  Van  der  Zlel"  has  extended 
the  solution  to  the  one-dimensional  p-n-p 
transistor  structure.  In  both  studies  the 
effects  of  surfaces  were  considered  in  an  ap¬ 
proximate  manner. 

>  However,  surface  conditions  have  been 
found  to  Influence  markedly  the  performance  of 
p-n  junction  diodes  and  transistors  (Kingston") 
Considerable  theoretical  work  has  been  done  to 
understand  the  signal  properties  (voltage,  cur¬ 
rent,  frequency  relations)  of  p-n  junction 
devices,  considered  as  three-dimensional  struct 
ures  (Shockley7,  Van  Roosbroeck") .  It  Is  the 
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purpose  of  this  thesis  to  develop  a  theory  of 
noise  which  considers  the  p-n  junction  as  a 
three  dimensional  system,  and  which  treats  the 
effects  of  surfaces  in  an  exact  manner. 

A  second  objective  of  the  thesis  is  to 
test  and  extend  a  powerful  method  developed  hy 
Petritz3>^  for  studying  complicated  random 
processes.  This  aspect  of  the  work  is  of  inter 
bst  in  the  general  theory  of  random  processes. 
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CHAPTER  II 

PROBLEM  AND  METHOD  OF  SOLUTION 


2.1  Introduction 

We  assume  that  semiconductor  noise  is  a  * 
stationary,  ergodic  and  Markovian  random  process. 
Considering  the  local  hole  density,  Pt(x,y,z,t) 
as  a  random  variable,  this  is  a  three-fold  infin¬ 
ite  random  process.  In  order  to  solve  such  a 
complicated  problem,  we  have  generalized  a  method 
used  originally  by  Petritz.*  This  method  employs 
the  Kolmogorov- -Fokker-Planck  (KFP)  and  the 
Langevin  techniques  to  describe  the  noise  (Fel¬ 
ler”,  Chandrasekhar10,  Uhlenbeck  and  Ornstein11, 
Wang  and  Uhlenbeck12). 


2.2  The  Kolmogorov-Fokker-Planck  Equation  Approach 


IP 

The  KFP  equations  for  the  three  dimensional 
semiconductor  are  given  by 
!£lTn0CrUTO(rl,t)  -  -  P(-m.t.T)|Tn(r),t)2  QMt>)|K M) 

KlW™  (1) 


*  2  P(m.W|KW,t)  Q(k(r)|-m(r))  . 

kjy  to 

P(nu/m,t)  1b  the  conditional  probability  of 
finding  the  random  variable  with  a  value  m  after 


the  time  t,  if  at  zero  time  the  random  variable 


had  a  value  m0.  The  random  variable  Is  the  hole 
density  In  the  n-type  semiconductor.  The  Bymbol 
r  represents  r(x,y,z),  a  function  of  the  three 
rectangular  coordinates.  Q  is  a  transition 
probability  and  is  defined^  by  the  equation, 


P(,Klim,&t)=  +  order  (A t)2- .  (2) 
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Q  describes  how  the  system  changes  In  an  Infin¬ 
itesimal  interval  of  time.  At,  and  characterizes 
the  stochastic  process.  For  the  semiconductor 
problem,  Q  is  independent  of  time  and  the  pro¬ 
cess  is  stationary.  Q  is  non-zero  and  less  than 
unity  and  the  process  is  ergodle. 

The  interpretation  of  equation  (l)  is  that 
the  rate  at  which  the  conditional  probability 
P(m0/m,t)  changes  with  time  results  from  transi¬ 
tions  away  from  and  to  the  desired  state.  Equa¬ 
tion  (1)  is  subject  to  the  boundary  condition 

P ,0) :s  &  -m,  ,-m  ,  (^) 

where  b,,,  wis  the  Kronecker  delta. 

Since  a  random  process  is  characterized  by 
transition  probabilities,  we  list  them  for  the 
p-n  junction: 


Q(m  lr)|  rn^-l)  =  m(r)//<  , 

bulk  recombination; 

(4) 

^(rn(r)|'m(r)4- 1)  - ^j)/^  , 

bulk  emission; 

(5) 

bulk  diffusion 
decrease ; 

(6) 

bulk  diffusion 
Increase . 

(7) 

At  the  transverse  surfaces,  the  surface  transi¬ 
tion  probabilities  are 

mOrO/tj  , 

surface  recombin¬ 
ation; 

(8) 

surface  emission. 

(9) 

rs  denotes  that  the  random  variable  Is  evaluated 
only  at  the  surface;  t  is  the  bulk  recombination 
lifetime  of  a  hole  in  an  excited  state;  td  and 
t3  designate  the  lifetime  of  bulk  diffusion 
and  surface  recombination  respectively; 
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and  <)  designates  ensemble  average.  For  a 
stationary  random  process,  time  and  ensemble 
averages  are  equal. 

The  KFP  equations  with  the  set  of  transi¬ 
tion  probabilities  given  in  equations  (l)  and 
(4)  to  (9)  comprise  a  three  fold  infinity  of 
differential  equations,  since  p(r,i)  depends 
continuously  on  x,  y  and  z.  We  have  not  solved 
this  complete  set  of  equations,  but  we  later  use 
some  KFP  equations  to  solve  the  noise  in  an 
infinitesimal  region  of  the  semiconductor. 


2.3  The  Langevin  Equation  Approach 


The  Langevin  equation  is  a  deterministic 
equation  of  a  system  excited  by  random  noise 
sources. For  a  particle  in  a  viscous  medium 
it  is 


!f  +  0u=Mt),  (io) 

where  u  is  the  velocity  of  the  particle,  is 
the  viscosity,  and  A(t)  is  the  random  force. 

The  two  assumptions  made  are  that  A(t)  is  inde¬ 
pendent  of  u  and  that  A(t)  varies  extremely 
rapidly  compared  to  the  variation  of  u. 

By  generalizing  the  above  concept  we  have 
a  suitable  method  for  solving  the  p-n  junction 
noise  problem.  The  deterministic  equations  for 
minority  carrier  flow  are; 7 >3 

SRfrQ  v  .  !_  v.T(r,n  =  o  , 

»t  ^  q.  J  (xi) 


jCr,t)=  ECr,i)“lDvpttnt)^-f>Ppt^(12) 


^fte  subscript  t  denotes  the,  total  hole  density, 

3  is  the  current  density, E  is  the  electric 
field  intensity,  is  the  hole  concentration 

at  at  time  is  the  hole  mobility,  p„ 

is  the  hole  concentration  at  thermal  equilibrium, 
t  is  the  mean  lifetime  of  a  hole  in  bulk  n-type 
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semiconductors,  D  is  the  diffusion  constant,  and 
q  is  the  electronic  charge.  We  assume  that  the 
diffusion  current  is  much  greater  than  the  con¬ 
duction  current.  By  introducing  appropriate 
noise  sources  into  the  above  equations  we  have 
the  three  dimensional  generalization  of  the 
simple  Langevin  equation  (10): 

ayfrvt)  +  .£.w»t)  +  Jr-v-fat)  =  Sr(v,t)  , 

d-t  «  ‘l  (13) 


Ja\t)  +  cjD  v  p  = SB  CP,t) .  (!4) 

S  and  SD  are  noise  sources.  The  variable  p  is 
the  deviation  of  the  hole  density  from  its  equil¬ 
ibrium  value,  like  Eq,  (Bl).  (A  letter  preceding 
an  equation  number  indicates  the  appendix  in 
which  the  equation  is  found.) 

It  is  important  to  note  that  there  exists 
no  a  priori  knowledge  of  the  noise  sources; 
their  solution  is  a  key  part  of  this  work. 

After  finding  expressions  for  these  noise 
sources,  we  solve  the  deterministic  Eqs.  (13)  and 
(14)  and  find  the  noise  spectrum  of  the  p-n  junc¬ 
tion.  The  latter  step  involves  the  use  of  scalar 
and  tensor  Green's  functions. 
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CHAPTER  III 

THREE-DIMENSIONAL  NOISE  SOURCES 
FOR  A  p~n  JUNCTION 


5,1  Introduction 

It  is  necessary  to  derive  explicit 
expressions  for  the  noise  sources  appearing  in 
equations  (13)  and  (14).  We  consider  the  anal¬ 
ogous  but  considerably  simpler  problem  of  a 
one-dimensional  transmission  line. 


5.2  The  Inhomogeneous  Transmission 
Line  Equations1^ 

The  homogeneous  differential  equations 
for  a  one -dimensional  transmission  line  with  no 
series  Inductance  are 

^t'---Ca^)-GV(x,«,  (15) 

--  -  RI  (x,t)  •  (l6) 

c3% 

I  is  the  current  flowing  in  the  line;  V  is  the 
voltage  across  the  line;  0,  0  and  R  are  the 
capacitance,  shunt  conductance  and  resistance 
per  unit  length  of  line  respectively.  An  infin¬ 
itesimal  section  of  line  is  shown  in  Figure  1, 
page  8.  Taking  the  Fourier  transform14  of  the 
voltage  and  current  (either  represented  by  F)  we 
have 

F(t)=  /  F(f)  exp(*wt)  ^ 

'oo  • 
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-  (iuoC  +  (jr)V  -  -YV  , 

(18) 

(19) 

Y  -  iu°  C  +  G  1 

(20) 

A 

uo-a'rt t  . 

(21) 

I  and  V  are  functions  of  frequency,  f,  and 
position,  x. 

We  now  consider  voltage  and  current 
sources  (generators)  in  the  line;  the  total 
source  in  the  series  arm  is  AV0  (xi)/Ax,  while 
in  parallel  with  the  admittance  Y  ’is  a  current 
source  AI0(xf)/Ax.  in  the  limit  as  Ax  ap¬ 
proaches  zerfo ,  we  have 


AX  -  ax 

(22) 

where  F0 

Io(*,f). 

mission 

(x,.f)  represents  either  VQ(x,f)  or 

The  resulting  inhomogeneous  trans¬ 
line  equations  are 

±L  +  ^v=  -  > 

ax  1  *  ax 

(23) 

4  V  ,  -n  ]  ~  _  §¥• 

OX  +  m  dX 

(24) 

3-3  Analogy  between  the  Transmission  Line 
and  the  One  Dimensional  p-n  Junction 

The  three  dimensional  p-n  junction,  Eqs. 
(13)  and  (l4),  reduce  to  one  dimension  when  the 
functions  considered  are  constant  in  the  y  and 
z  directions.  Using;  Fourier  transforms,  Eqs. 
(lj)  and  (l4)  become 
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-Ik  • 

(25) 

038.  +  4-  -  “  ’ 

a* 

(26) 

(1+  io)<)/D  ^  • 

(27) 

Here  j  and  p  are  functions  of  frequency  and  x. 
The  quantities  in  Eqs.  (25) and  (26)  analogous 
to  those  in  Eqs.  (23) and  (24)  are 


I~j  ;  ^P  •>  (28) 

Y-v  q^OJ  +  <'');  R~  \/*iP  >  (29) 

al./a*^  d>/a7L  >  (50) 

aV0/ax  sp./9^  •  (51) 


5.4  The  Three  Dimensional  Inhomogeneous 
Differential  Equations  for  a 
Semiconductor  with 
a  p-n  Junction 

Since  the  current  density  Is  a  vector  and 
the  excess  hole  density  is  a  scalar,  there  are 
four  more  equations  like  (25)  and  (26)  for  the 
y  and  z  directions.  These  six  equations  can  be 
written  as 


V-^  +  p  -  > 

(52) 

V7p  +  3/^  =  ~^P° 

(55) 

10 
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When  the  variables  are  separated,  we  get 


72p-KV  -  *5  ’ 

(34) 

s=3>S.=  ^p.-5-j../4d, 

(35) 

(36) 

+3*'  vp0  . 

(37) 

5.5  Discussion  of  the  Random  Noise  Sources 

In  the  first  order  differential  equations 
(32)  and  (33)  there  are  two  sources  of  noise: 
v-jn  and  vp.  .  The  first  is  the  divergence  of 
the  hole  current  and  is  a  scalar  source,  This 
is  due  to  hole  recombination  with  electrons. 

The  second  noise  source  is  the  gradient  of  the 
hole  density,  vp,  ,  and  is  a  vector  source.  T^he 
diffusion  current  density  is  proportional  to  vp0, 
and  is  in  the  direction  from  greater  to  lesser 
hole  density. 

In  addition  to  the  above  noise  sources 
there  are  sources  which  result  from  the  recom¬ 
bination  or  emission  of  holes  in  surface  states. 
This  noise  is  related  to  the  flow  of  holes  into 
the  surface  and  is  directed  normal  to  the  sur¬ 
face. 

The  modified  noise  sources,  S^,  and  , 

Eqs.  (35)  and  (37);  are  not  new  sources  but 
result  from  mathematical  operations  on  the  physi¬ 
cal  sources. 


3.6  Method  of  Deriving  the  Noise  Sources 

We  derive  explicit  expressions  for  the 
noise  sources  following  the  method  first  used  by 
Petritz.^’^This  method  uses  the  KFP  equation  to 
solve  for  the  spectrum  of  the  noise  in  an 
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.Infinitesimal  region  of  the  semiconductor . 

Then  an  appropriate  deterministic  equation 
with  unknown  noise  sources  is  set  up  for  the 
same  infinitesimal  region  (local  Langevin 
equation).  Knowing  the  spectrum  from  the  KFP 
solution,  one  is  able  to  derive  expressions  for 
the  noise  sources  appearing  in  the  local  Lange¬ 
vin  equations.  These  sources  turn  out  to  be 
appropriate  for  use  in  the  Langevin  equations 
of  the  whole  semiconductor,  Eqs.  (24)  to  (27). 


2-7  Models  Used  for  Determination 
of  Noise  Sources 

The  determination  of  the  noise  sources 
is  simplified  because  of  the  assumption  of 
statistical  independence  of  the  various  elemen¬ 
tary  processes.  To  set  up  a  model  which  isolates 
each  source,  we  cut  the  three  dimensional  semi¬ 
conductor  into  infinitesimal  cubes  without 
changing  their  hole  density,  and  apply  the  appro¬ 
priate  boundary  conditions.  For  recombination 
the  cube  is  an  Interior  one  with  perfectly 
reflecting  boundaries.  The  charge  density 
remains  uniform  throughout;  the  only  decay  is 
due  to  the  bulk  recombination  time  constant 
since  diffusion  currents  require  gradients. 

For  the  bulk  diffusion  sources  the  cube 
is  an  interior  one  with  perfectly  absorbing 
boundaries  on  two  opposite  faces  and  perfectly 
reflecting  boundaries  on  the  other  faces.  The 
volume  is  so  small  that  the  concentration  grad¬ 
ients  cause  large  diffusion  currents  while 
relatively  few  holes  are  lost  by  recombination. 

For  the  surface  recombination  noise  source, 
the  cube  is  at  the  surface  of  the  semiconductor. 
It3  dimension  perpendicular  to  the  surface  is 
very  small  compared  to  the  others  and  the  bound¬ 
aries  are  perfectly  reflecting  except  for  the 
original  semiconductor  boundary.  In  this  sur¬ 
face  element  the  dominant  process  is  recombina¬ 
tion  and  emission  from  the  surface  states. 
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Because  of  the  thinness  of  the  element,  bulk 
recombination  and  diffusion  are  relatively 
unimportant. 


5.8  The  Bulk  Recombination  Noise  Source 

The  excess  hole  density  Langevln  equation 
for  bulk  recombination  and  emission  is  obtained 
from  the  model  in  Section  3.7  and  from  Eqs.  (10) 
and  (13), 

ft  +  f ;  ^  •  (38) 

The  random  force  is  frj/q  .  The  time  constant 
for  bulk  renomh-1  nation  r  determines  the  transi¬ 
tion  probabilities j 

Q>(N|N“I)=N/'£  ,  bulk  hole-electron  (39) 

recombination; 

Qr  (N |N+ 1)  - ,  excitation  of  a  (40) 

hole . 

Here  N  and  <N>  are  the  total  and  the  average 
number  of  holes  in  the  cube  respectively.  The 
minority  charge  density  is  assumed  uniform  over 
the  volume  v  of  the  cube,  thus 


S  s  V  pt  .  (41) 

Substituting  these  transition  probabilities  into 
the  KFP  equation  (l),  we  get  for  the  small  vol¬ 
ume  under  consideration, 

.ipwu-psqw  +p(R|s-i,t)  m _  („a) 

9"t  lX*  X  '  ' 

+  P(W.|N+l,t)  [M^J]  . 

Multiplying  equation  (42)  by  N,  the  number  of 
holes  at  time  t  if  NQ  existed  at  the  initial 
time,  and  summing  over  the  ensemble,  there 
results 
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a  <Nt>.  /at  =  -  <Ni>y  -t  ,  (43 ) 

N^N-<N)  ,  (44) 

and  where  the  conditional  average  for  N  is 
defined  as 

<N>.  =  J  N  .  (“5) 

The  solution  of  equation  (43)  is 

<K>o  =  Nie  expC-t/O  ,  (46) 

where  N-^q  is  defined  with  the  aid  of  Eq.  (3)  as 

<Nj)0=  Ni0  »  t*0  .  (47) 


This  conditional  average  is  transformed 
into  the  correlation  function  using^ 


(48) 


with  equation  (45)  this  expression  becomes 

p(tl=j]N„  W*U  <N,>.  .  (49) 

*10 

w(Nlo)  is  the  probability  of  Nj.o»  When  Eq.  (46) 
is  substituted  into  (49),  the  correlation  func¬ 
tion  becomes 

p(f)»  (N’,o)’  exp("V^)  •  (50) 

The  Wiener-Khintchine  Theorem-^  transforms 
the  correlation  function  into  the  spectrum 
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w(£)=-4{  pft)  cos  cot  dt  .  (51) 

When  equation  (50)  is  substituted  into  equation 
(51)  and  the  result  integrated,  there  results 


(52) 

K2  is  defined  by  Eq.  (27)  and  the  symbol  ||  means 
absolute  value. 

Equation  (52)  is  the  spectrum  of  the  total 
hole  variation.  However,  we  desire  an  expres¬ 
sion  for  the  noise  source  which  appears  in  the 
Langevln  frequency  equation  (54).  We  can  solve 
for  the  spectrum  of  the  source  now  that  we  have 
determined  the  spectrum  of  the  total  hole  varia¬ 
tion.  This  is  the  inverse  of  what  is  normally 
done  with  the  Langevln  equation.  The  normal 
procedure  is  to  postulate  a  white  noise  source 
and  to  solve  for  the  fluctuations  in  the  total 
hole  density.  The  magnitude  of  the  white  noise 
source  is  determined  by  considerations  of  statis¬ 
tical  mechanics  and  thermal  equilibrium.  Since 
we  are  interested  in  non-equilibrium  as  well  as 
thermal  equilibrium  noise  sources,  we  cannot 
use  this  normal  procedure.  Instead  the  local 
KFP  equation  is  used  to  determine  the  thermal 
equilibrium  and  nonequilibrium  spectrum  of  the 
total  hole  variation,  and  the  noise  sources 
derived  have  general  validity.  This  combined 
use  of  the  KFP  and  the  Langevln  methods  locally 
was  first  done  by  Petrltz*  and  appears  to  be  a 
powerful  technique  for  solving  complicated  ran¬ 
dom  processes. 

We  use  this  technique  to  transform  the 
spectrum  of  the  total  hole.1  variation,  Eq.  (52), 
into  the  noise  source  associated  with  the  recom¬ 
bination  dissipative  process.  This  is  done  by 
integrating  equation  (58)  over  the  volume  of  the 
cube;  inserting 


15 


HAVORD  Report  5762 


N^JpJrtdv  ,  (53) 

where  p-^r)  is  like  p.^0),  Eq.  (Bjj);  using 

Ir=  J  (?■!')  5  (54) 

and  using  Eq.  (17).  This  yields 

I„=  qCK'-N.  .  (55) 

The  spectrum  of  Ir  is 

^(IT/)=  ^IK4!  w(IN,f).  (56) 

Equation  (52)  substituted  into  Eq.  (56)  obtains 

.  (57) 


This  is  the  noise  source  associated  with  the 
recombination  dissipative  process.  It  is  inde¬ 
pendent  of  frequency  and  therefore  1b  a  white 
noise  source. 

Equation  (57)  Is  not  the  recombination 
noise  source  in  equation  (55).  This  source  is 

wos/)=v*r(i|ifr).  (58) 


A  relation  between  Eqs.  (57)  and  (58)  is 
determined  with  the  equation 


Lim. 

ftV-*  0 


and  the  generalized  Impedance  theorem. 


(59) 
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The  impedance  theorem  (Lawson  and  Uhlen- 
beck1-})  states  that  by  multiplying  the  current 
noise  spectrum  Wj_  in  a  linear  system  by  the 
absolute  value  squared  of  its  impedance  function, 
Z(f)Z*(f),  wj_  is  transformed  into  the  voltage 
spectrum  wc.  In  symbols 

wo=  ZZ*.  (60) 

We  generalize  this  impedance  theorem  by 
letting  Wj  and  wQ  be  any  two  spectra  in  a  linear 
system  which  are  related  by  a  factor  of  propor¬ 
tionality  Z(f)Z*(f). 

The  explicit  expression  for  the  recombin¬ 
ation  source  spectrum  is 


w(|Srl*)=  4  <<)Ad’  4v\  (61) 

To  express  the  recombination  source,  Eq. 
(61),  in  known  parameters,  <N*>  is  evaluated 
by  multiplying  Eq.  (42)  by  N*  and  summing;  the 
result  is 

*  Ci W.  *  <N) -  +  <n>.)A . . 

at  (62) 


When  the  time  of  observation  of  N  goes  to  infin¬ 
ity, 


<*’>  =  (N1)-  <*»>  (6?) 

since 

<N>,  =  <N>-  (64) 

With  Eqs,  (65)  and  (41)  Eq.  (61)  becomes 

vr (iS^l*)  =  4  <pt(.r)>/^  D*  dv),  (65) 
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where  <p^(r))  is  given  by  Eq.  (Bll).  This  16 
the  desired  noise  source  for  bulk  recombination. 


5.9  The  Diffusion  Noise  Source 

Solving  Eqs.  (13)  and  (14)  for  p, 
letting  Sp  be  zero  and  using  the  model  for  the 
diffusion  noise  source.  Section  (3.7),  we  get 

5v  '  (66) 


where  x^  is  x,  y,  or  z.  We  assume  that  diffusion 
in  the  three  directions  is  statistically  inde¬ 
pendent.  Equation  (66)  becomes  a  set  of  three: 


where  u  stands  for  x,  y  or  z. 

To  write  equation  (66)  in  the  Langevin 
form,  the  spacial  term  is  transformed  to  contain 
a  time  constant.  We  write  the  second  derivative 
for  the  finite  but  small  cube  and  use  the  densi¬ 
ties 


p(0)=0  ,  pCa-a^ 0  , (6Q) 


In  this  differentiation  the  diffusing  direction 
is  u  and  the  length  of  the  cube  in  this  direction 
is  Au.  The  second  derivative  is 


a_E  -  _  JLp 
au*  iu1 


and  equation  (67)  becomes 


(69) 
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where  Auy2D  .  (71) 


With  these  time  constants  the  transition 
probabilities  for  the  u  direction  are 


,  loss  by  diffusion;  (72) 

Q»tt0$lNp+l)  ,  gain  by  diffusion.  (73) 


Here  Np  is  the  total  number  of  holes  in  the  cube 
with  dimensions  Ax,  Ay,  Az  and  with  the  diffu¬ 
sion  boundary  conditions: 


=  Jp  dv  =  AX  ^  AT.  p/2.  .  (74) 

Using  the  techniques  of  the  previous  sec¬ 
tion,  the  noise  source  for  diffusion  in  the  u 
direction  is 


wOS,(«\h-  lK4KcfO<p.(rt>/«  ^  a*  .  (75) 


3,10  The  Surface  Recombination  Noise  Source 

With  the  model  for  the  surface  recombina¬ 
tion  noise  Bource,  Section  3.7,  the  Langevln 
equation  (10)  becomes 

+  IgjL  a  Sil.  .  (76) 

The  pia  is  evaluated  at  the  semiconductor 
boundary  and  the  is  the  current  flowing  in 
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th r>  direction  normal  to  the  surface  i.  These 
surface  recombination  time  constants  define  the 
transition  probabilities: 

by  recombination  at  the  surface  ; 

by  emission  at  the  surface  . 

Ns  is  the  total  number  of  minority  carriers 
evaluated  at  the  surface  if  the  surface  layer 
were  Axi  thick: 


(77) 

(78) 


1*3=  /pis'  ATCv-pis  •  (79) 


Using  the  techniques  demonstrated  in  Sec¬ 
tion  J.7j  the  surface  recombination  noise  source 
is 


vte  01%/)  =  4  <?t»>  a /tf  dx  ,  (80) 

where  1  and  k  are  either  the  pair  y,  z  or  z,  yj 
s  is  the  surface  recombination  velocity  and  is 
assumed  to  be  the  same  on  the  y  and  z  surfaces. 
The  relation  for  the  surface  recombination  vel¬ 
ocity  is  given  by  Rittner,-*-© 


Si= i  .  (81) 

The  rate  of  surface  recombination  acts  as  if  a 
current  of  holes  were  drifting  into  the  surface 
with  an  average  velocity  s  and  being  removed. 1 


3.11  Noise  Source  and  p-n  Junction  Equations 
The  two  bulk  sources  in  equations  (65)  and 
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(75)  are  expressed  per  unit  volume.  Since  the 
differential  equations  are  for  a  small  volume 
of  the  three  dimensional  semiconductor,  the  noise 
sources  are  in  the  correct  form.  That  these 
sources  correspond  to  physical  processes  was  also 
shown  by  analogy  with  the  transmission  line. 
Section  3.3. 

When  the  variables  are  separated  in  the  two 
first  order  semiconductor  equations,  spacial 
differential  operations  are  performed  on  some 
of  the  noise  sources.  This  changes  their  nature 
from  those  calculated  in  Sections  3.8  and  3.9. 

In  the  equations  for  excess  hole  density  (34) 
and  (35)  the  nature  of  the  recombination  source 
Sr  is  not  changed  and  these  equations  can  be 
solved  for  the  recombination  noise  spectrum. 
However,  the  diffusion  source  SD  has  been  dif¬ 
ferentiated  spacially  and  becomes  Sr,.  By  using 
the  vector  analogue  of  integration  by  parts,  the 
diffusion  noise  souroe  transforms  to  the  correct 
form.  This  technique  is  used  in  Chapter  VII. 

In  the  diffusion  current  density  equations 
(3 6}  and  (37)  the  diffusion  source  is  still  Sn. 
These  equations  can  be  solved  for  the  diffusion 
noise  spectrum  with  the  tensor  Green' 3  function. 
This  approach  is  followed  in  Chapter  VI. 
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CHAPTER  IV 

SCALAR  INHOMOGENEOUS  SEMICONDUCTOR  EQUATION 
AND  GREEN'S  FUNCTION 


4.1  Formal  Solution  of  the  Scalar 
Inhomogeneous  Equation 

Having  derived  explicit  expressions  for 
the  noise  sources,  we  consider  the  Inhomogeneous 
differential  equations  for  the  semiconductor  with 
a  p-n  Junction,  Eqs.  (34)  to  (37).  Green's  func¬ 
tions  are  useful  for  solving  inhomogeneous  partial 
differential  equations.  In  our  problem  we  use 
both  scalar  and  tensor  Green's  functions,  the 
latter  because  of  the  vector  nature  of  the 
sources  in  Eq.  (36).  A  summary  of  the  proper¬ 
ties  of  scalar  and  tensor  Green's  functions  is 
given  in  Appendix  A. 

The  formal  solution  of  Eq.  (34)  in  terms 
of  a  scalar  Green's  function  is,1'* 


■pW * J [GftpW) - P  W ^ ^  +  P ^  ^  •  (82) 


Throughout  the  paper  the  zero  subscript  denotes 
the  source  coordinates,  while  the  coordinates 
without  subscripts  are  the  observation  ones. 

The  surface  integral  gives  the  contribution  for 
noise  sources  at  the  surfaces,  while  the  volume 
integral  is  for  the  volume  sources. 
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4.2  General  Discussion  of 
Scalar  Eigenfunctions 

We  construct  an  explicit  expression  for 
the  scalar  Green's  function  in  terms  of  a  series 
of  scalar  eigenfunctions.  An  eigenfunction  is 
the  solution  of  an  ordinary  homogeneous  differ¬ 
ential  equation  containing  a  separation  constant 
which  satisfies  simple  boundary  conditions.  The 
values  of  the  separation  constants  which  allow 
the  eigenfunction  to  fit  the  conditions  are 
called  eigenvalues.  In  physics  it  is  assumed 
that  the  Dirac -delta  function  related  to  our 
Green's  function  can  be  expanded  in  terms  of  a 
complete,  orthogonal  set  of  eigenfunctions. 

The  orthogonality  condition  for  a  set  of 
eigenfunctions  is 


^  —  ^.-wvn  •  (8?) 


Here  is  the  Kronecker  delta  and  Llm_ 

is  the  normalization  constant.  The  differential 
equation  the  eigenfunctions  of  this  problem  must 
satisfy  is 

0  •  (84) 


2. 

Klmn  13  the  separation  constant  and  specific 
values  of  the  separation  constant  for  which  the 
above  equation  can  be  solved  are  the  eigenvalues. 


4.5  The  Eigenfunction  Expansion  of  the  Scalar 
Green's  Function  for  Arbitrary  Surface 
Recombination  Velocity 


The  scalar  Green's  function  is  now  expanded 
in  a  series  of  scalar  eigenfunctions: 


C--Y, 


n 


(85) 
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This  series  Is  substituted  into  equation  (Aljj), 
noting  that  the  coefficients  are  not  factions 
of  r.  Multlplyijig  by  Fr^Cr).  integrating  over 
the  volume,  and  using  Eqs.  (83)  and  (84),  we  find 

/A»Jlw.Tl.^K*TrvX>.  (86) 

and  ,  (8?) 

where  K*+  vi if  *  (88) 


This  Green1 s  function  satisfies  the  reciprocity 
condition  since  it  is  symmetrical  in  the  source 
and  observation  coordinates. 

The  semiconductor  geometry  is  shown  in 
Figure  2.  The  p-n  junction  is  located  at  the 
x=0  plane  and  the  origin  of  coordinates  lies 
at  the  center  of  this  face.  The  rectangular 
parallelepiped  is  bounded  by  the  plane b  x&O, 
x=a,  y=b,  y=-b,  z=c,  z*-c. 

The  three  dimensional  rectangular  coordin¬ 
ate  system  is  a  separable  system  and  the  eigen¬ 
function  can  be  written  as  the  product  of  three 
factors: 


Each  factor  is  the  eigenfunction  which  satisfies 
the  boundary  condition  in  one  coordinate.  Fur¬ 
thermore  Eq.  (84)  separates  into  three  equations: 

,  (90) 


where  r  is  any  index,  1,  m  or  n. 

The  boundary  conditions  in  the  x-dlrectlon 
are  that  p*o  at  x=o  and  x=a.  This  implies  an 
ohmic  contact  at  x=a  and  short-circuited  condi¬ 
tions  at  the  x=0  and  x=a  planes.  The  eigenfunction 
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Fig.  2.  The  geometry  and  coordinate  system  for 
the  three  dimensional  semiconductor  with  finite 
surface  recombination  velocity  on  the  transverse 
surfaces  y=b,  y=-b,  z=c,  z=-c.  The  p-n  Junction 
is  at  the  x=0  boundary  and  an  ohmic  contact  Is 
at  the  x=a  boundary. 
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which  satisfies  these  conditions  is 


Fj,-  si*  rtW3-  •  (91) 

On  the  y  and  z  boundaries  current  moving 
into  the  surface  is  proportional  to  the  excess 
hole  density.!  The  constant  of  proportionality 
is  the  surface  recombination  velocity  s: 

:»/<!=  tsP>  .  (92) 

The  coordinate  u  stands  for  either  y  or  z  and  * 
stands  for  the  y  or  z  boundary  surfaces.  Using 
the  homogeneous  form  of  Eq.  (14),  Eq.  (92) 
becomes 

ap/aui=?  Sp/D  ,  (95) 

Equations  (90)  and  (93)  have  the  following 
solutions:  In  the  y-directlon  the  cosine  and 
sine  eigenfunctions  are 


qo*  i  IjJ  -  sift  Pm y*  (9^0 

In  the  z-direction  the  cosine  and  sine  eigen¬ 
functions  are 


Fw»  co$  ,Fy  =  fljn  ’  (95) 

The  boundary  condition,  Eq,  (93),  becomes 
for  the  cosine  eigenfunctions: 

f*  *  'V®  »  (96) 


and  for  sine  eigenfunctions: 
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fif  cot  JP  =  -  3 Vc  •  (97) 

The  subscript  r  stands  for  either  m  or  n,  the 
subscript  p  stands  for  either  or  \>  ,  and 
stands  for  either  *  b  ono,  respectively.  (The 
word  "respectively"  in  this  expression  means 
that  only  the  values  (m,  M  ,  t. b)  or(n,y  ,  *c) 
can  occur  together  in  the  above  equations^  We 
assume  that  the  surface  recombination  velocity 
is  the  same  for  an  opposite  pair  of  surfaces. 

Combining  Eqs.  (91)*  (9^+)»  and  (95)  the 
eigenfunctions  for  the  rectangular  parallelepiped 
semiconductor  are 

Pjnul,  **  ;x(/W)  SiiCA.*).  os) 

where  all  possible  sine  and  cosine  combinations 
are  taken.  The  symbols  M  and  N  stand  for  a  y 
or  z  index  respectively.  These  eigenfunctions 
are  complete  since  they  satisfy  Eq.  (84)  and 
the  homogeneous  boundary  conditions.14  The 
scalar  Green's  function  Eq.  (87)  satisfies  the 
same  boundary  conditions  as  the  eigenfunctions. 
The  eigenfunctions  are  orthogonal;  from  Eqs. 

(83)  and  (98)  we  find 

dv  s  (99) 

and 

(100) 

while  any  mixing  of  the  m,p  and  n ,i>  give  a 
null  integral.  Since  the  eigenfunctions  are 
separable  we  can  integrate  the  various  coor¬ 
dinate  integrals  separately; 

(101) 
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J*  ,«jja  Sin  dK  ~  &  ^>,4'  •  (102) 

/%o?  £„u  co?  j^.-u  du  =  «l'+  )^t,r~  L*r  6*\r',(i02) 

(*n/S,u  ™^-U  du-O-^IV^/V.  (104) 


f*- cat  fit*  5in^udu=0.  (105) 


The  r  stands  for  m  or  n  While  p  stands  for_>^ 
or  )?  . 

Substituting  Eq.  (98)  into  (84),  the 
eigenvalues  are 


KUut*Ki -'•*0  tin  ~  0r)  +£  'V3*  ‘ 


(106) 


If  we  let  a,  the  x-dimenslon  of  the  para¬ 
llelepiped,  go  to  infinity,  the  discrete  sum 
goes  over  to  an  Integral  and  can  be  integrated. 
To  carry  out  this  limit  operation,  Green's  func¬ 
tion  is  written  as 


rr..WF.>l  VI 
U  L«K 


„■  -nk  .uKSXo 
Sir  -g—  an  -a-* 

S*  -t-Kttwn 


•  (107) 


The  variables  H  and  ATJ  take  on  the  values  ^Vd 
and  fl/a,  respectively,  and  are  put  into  Eq.  (107) 
When  a  goes  to  infinity,  Eq.  (107)  becomes 


(n)Wr) 


KtW 


JJrtK 


/ ; 


din  rx  sin  d* 

KlK«»  ’  (108) 
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When  we  integrate  equation  (108),  Green's 
function  becomes 


1  +  i£L 

*  D 


(112) 


is  defined  by  Eqs.  (101).  (10J)  and  (104) 
ana  by  Eqs.  (106)  and  (88).  The  plus  signB 

belong  to  the  m,n  indices,  while  the  minus  signs 


to  the^,*?  Indices.  This  is  the  desired  form 
for  the  scalar  Green's  function  for  a  rectangular 


parallelepiped  semiconductor  when  the  x-dimension 
goes  to  Infinity.  The  assumption  that  a  goes  to 
infinity  means  physically  that  the  diffusion 
length  of  the  minority  carriers  is  much  less  than 
a.  Therefore  minority  carriers  injected  at  x=o 
will  recombine  before  they  reach  the  contact  at 
x=a.  This  assumption  is  valid  for  p-n  Junction 
diodes.  In  the  case  of  transistors,  a  second 
p-n  junction  is  at  the  distance  a,  and  in  this 
case,  the  diffusion  length  is  much  greater  than 


a. 


Our  solution  is  therefore  directly  appli¬ 
cable  to  p-n  junction  dlodeaj  transistors  can  be 


29 


NAVORD  Report  5762 

handled  in  the  same  general  way,  but  a  is  kept 
finite  and  the  boundary  condition  at  a  is 
changed.  For  simplicity  we  confine  this  work 
to  the  p-n  junction  diode  and  assume  a  is 
infinite. 


4.4  The  Eigenfunction  Expansion  for  the 
Scalar  Green's  Function  for  Infinite 
Surface  Recombination  Velocity 


Green's  function  for  the  case  of  infinite 
surface  recombination  velocity  on  the  transverse 
surfaces  and  zero  hole  density  on  the  longitud¬ 
inal  surfaces  is  called  the  infinite  s  case.  It 
is  convenient  to  change  the  coordinate  system 
from  that  of  Figure  2  to  that  of  Figure  in 
which  the  origin  is  at  one  corner  of  the  p-n 
junction  faDe  rather  than  at  its  center.  To 
designate  that  a  symbol  pertains  to  the  case  of 
infinite  s,  we  affix  the  superscript  «•>  to  the 
symbol . 

On  the  Jxmgitudinal  surfaces  the  boundary 
conditions  are  not  changed.  On  the  transverse 
surfaces  s  is  infinite.  From  Eq.  (93),  we  get 


JL  -  +  -P- 

s  “  +  n«t 

■a  ix. 


U*  -  <5 


(113) 


Here  <r  1b  the  u-boundary  surface:  0,  B  for  y 
or  0,  G  for  z.  With  s = »  ,  either  p  is  zero 
or  Sp/au  is  infinite.  The  last  relation 
requires  an  infinite  surface  current  and  is 
physically  impossible.  Therefore  p  is  zero  on 
the  surfaces  which  have  infinite  surface  recom¬ 
bination  velocity.  Thus  the  excess  hole  density 
is  zero  on  all  surfaces. 

Following  the  same  procedure  as  for  the 
case  of  arbitrary  s,  we  find 

<("=££  LZzMLsmW  [«f(-C.\»-*.V .Mj}  (n4) 

™  k: _ 
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Fig.  3-  The  geometry  and  coordinate  system  for 
the  three  dimensional  semiconductor  with  infin¬ 
ite  surface  recombination  velocity  on  the  trans¬ 
verse  surfaces  y=0,  y=B,  z=0,  z=C.  The  p-n 
Junction  is  on  the  x=0  boundary  and  an  ohmic 
contact  is  on  the  x=a  boundary. 
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r  «o  •  n 

V  -  Sin  VI 


sin  TVp. 


2. 


£. 

z 


KL  -  "•->  ^  (3rf+  (r^)  ■ 


(115) 

(116) 
(117) 


Equation  (114)  Is  the  scalar  Green's  function 
for  infinite  surface  recombination  velocity  on 
the  transverse  surfaces. 
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CHAPTER  V 

TENSOR  GREEN'S  FUNCTION  FOR  THE  SEMICONDUCTOR 
WITH  A  p-n  JUNCTION 


5.1  Formal  Solution  of  the  Diffusion  Current 
Density  with  the  Tensor  Green's  Function 

In  deriving  the  noise  from  a  p-n  junction 
the  vector  Inhomogeneous  differential  equation, 
Eq.  (56)j  must  be  solved;  one  method  of  solution 
is  with  a  tensor  Green's  function.  The  general 
properties  of  tensor  Green's  fupction3  are 
reviewed  in  Appendix  A. 

We  solve  the  inhomogeneous  current  density 
equations  (36)  and  (37)  formally  with  the  tensor 
Green's  function,  P,  defined  by  Eq.  (A14).  Equa¬ 
tion  (36)  is  postmultlplied  with  T  while  equation 
(A14)  Is  premultiplied  with  J  and  the  two  result¬ 
ing  expressions  are  subtracted.  The  following 
tensor  Identities  are  used: 


.  (us) 

and  v[j(?nl=(yT)('-rM'F?'P,).  (119) 


The  resulting  equation  is  integrated  over  the 
volume  using 


•F=  F  •■Jl  -  F 


(120) 
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for  any  vector  function  F  and  for  the  idem- 
factor  <2,  Eq.  (A10).  The  symbols  r  and  rQ  are 
Interchanged  and  the  reciprocity  condition,  Eq. 
(A22)  is  used.  The  current  density  is 

•j  w  -  /£$,•  j i cr.^j  p  -  T(oK*r]}  •  jtsfcynU. .  (121) 


5.2  Formal  Vector  Eigenfunction  Expansion 
of  the  Tensor  Green* s  Function 

In  order  to  use  Eq.  (121)  we  must  find  an 
explicit  expression  for  the  tensor  Green's  func 
tion.  This  la  accomplished  formally  with  a 
complete,  orthogonal  series  of  vector  eigen¬ 
functions,  'j.  .  These  vector  eigenfunctions 

must  satisfy  the  equation 


and  the  orthogonality  condition 

/ dv  “A-JU*  .  (123) 

We  expand  the  tensor  Green’s  function  in 
a  series  of  these  vector  eigenfunctions, 

r-lLl-W'  (124) 

Jj-WVl 

This  series  is  put  into  Eq.  (A14)  and  both  sides 
f  the  resulting  expression  are  multiplied  by 
tWn'and  integrated  over  the  volume.  Using  Eqs, 
123)*  (120),  and  (88)  the  vector  coefficients 
are 

•  (125) 
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Substituting  Eq.  (125)  into  (124)  the  tensor 
Green's  function  becomes 

r Ja-v,  <r)  3Wir.)  /A-*wv^ksi —  .  (126) 

A“VY\V' 

This  is  the  formal  expression  for  Green's  ten¬ 
sor  in  terms  of  a  complete  orthogonal  set  of 
vector  eigenfunctions. 


14 

5,5  Vector  Eigenfunctions 

The  vector  eigenfunction  solutions  of  the 
vector  equation  (56)  are  obtained  from  scalar 
eigenfunction  solutions  of  the  corresponding 
scalar  eigenfunction  equation  (84).  These  vec¬ 
tor  eigenfunctions  are  written  as 

,  (127) 

where  L ■, 


and  X,  Y,  Z  or  y  are  eigenfunction  solu¬ 

tions  of  equation  (84);  is  the  vector  normal 
to  the  surface  and  w  is  a  function  of  the  coor¬ 
dinate  in  the  direction  a,  .  The  eigenvalues, 
Klmn>  and  the  relative  magnitudes  of  X,  Y,  Z  or 
,%  are  adjusted  to  satisfy  the  boundary 
conditions. 

In  the  p-n  Junction  the  diffusion  current 
is  proportional  to  the  gradient  of  the  hole 
density.  Therefore  a  vector  eigenfunction  is 
the  gradient  of  a  scalar  eigenfunction.  The 
curl  of  these  vector  eigenfunctions  is  zero,  and 
M  and  N,  Eqs.  (127)  and  (128),  are  zero. 

For  an  arbitrary  surface  recombination 
velocity  the  vector  eigenfunctions  are 
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>.J( 


in..— Ol 


(129) 


by  (106) 


where  Pimn  are  defined  by  Eq.  (98).  Klipn  by  ( 
and  the  boundary  conditions  by  (96)  and  (97) • 
These  eigenfunctions  satisfy  the  differential 
equation  (84). 


5.4  Non-orthogonality  of  the  Vector 
Eigenfunctions  for  Finite  s 

To  investigate  the  orthogonality  of  the 
vector  eigenfunctions  we  write  one  of  them 
in  component  form: 


(130) 


T*  ^ 

M*»  f>M  J 

P„i. 

s  in  fa  t 

sm  *0— 

<T>  SI"  P**3 

COS  fa  y 

cos  ^,1 

si*  Po  ’>■ 

firi 

COj  Pm  3 

»'*  Pm  y 

whe:  M  refers  to  m  or,A  and  N  to  n  or  V  .  All 
possxble  combinations  of  the  cosine  and  sine 
factors  are  implied  in  each  term.  A  component 
of  J*imn  denoted  by  the  coordinate  written  ao 
a  superscript. 

The  orthogonality  of  the  set  of  vector 
eigenfunctions  is  determined  from  Eq.  (123): 
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To  evaluate  these  Integrals, Eqs.  (102)  to  (105) 
and  the  following  expressions  are  used: 


j  sin  6rxjL  sin 

4*nC.(/$r-^r-V] 

Sir  P,* 

Si*  (?«.■"< 

(132) 

_  jjv  **' 

•  rtr1 

fsitlVu  siin^'U  du  -  **  [\~ 

(sir  = 

;  r-v>’ 

(133) 

L 

W  J 

A 

l‘°4U  C0S/fU 

'  cos$°<  cosj^.«< 

(134) 

=  \lfJ> 

i  f*/ 

(c.o^U 

cefc^.u  du  ~  *  (* 

+  Sin  ifrA 

'A«  } 

- 1  +  ' 
~L/'  )/7 

(135) 

^ru 

-  «* 

coajSpiL  -  0 

(136) 

The  superscript  plus  is  not  written  unless 
there  is  a  chance  of  ambiguity.  Here  r  stands 
for  m  or  n,/>  stands  for  /*■  or  J  ,  and  u  for  y 
or  z,  respectively. 

The  components  of  Eq.  (ljl)  are  now 


u 


* 


uv 


5, 


(l?7) 
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The  y  component  has  two  cases: 

Case  1:  ,  M*M«,  N=N» 

/)!„„  . 


(138) 


Case  2:  M«M>,  N°N ' 

/i»MW  Ivm'n'  Op  Lm  H'wW  •  (^29) 


The  z-component  also  has  two  cases: 
Case  1:  i=l' ,  M«M',  N*N‘ 


/&.  &*,  ■»«  =  0><\  f~)\  0,  V,w  • (wo) 


Case  2:  L-l' ,  M«M*,  N»N« 


M  refers  to  either  m  or  and  N  refers  to  n  or 
v  .  The  5i,  j  are  Kronecker  deltas.  When 
two  symbols  are  written  In  a  column,  the  upper 
symbol  refers  to  either  m  or  n  and  the  lower  one 
to  m  or  v  .  The  Integrals  In  Eqs.  (138)  and 
(140)  do  not  vanish  even  when  M*M'  or  N^N';  the 
set  of  vector  eigenfunctions  Is  not  orthogonal. 
Therefore  we  are  not  able  to  construct  a  tensor 
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Green's  function  from  these  vector  eigenfunc¬ 
tions. 


5.5  Vector  Eigenfunctions  for  Infinite 
Surface  Recombination  Velocity 

In  the  special  case  of  infinite  surface 
recombination  velocity  a  complete  set  of  ortho¬ 
gonal  vector  eigenfunctions  can  be  derived  for 
the  three  dimensional  p-n  .junction.  These' eigen- 
.functionfc  are  fouhd  from  the  products  of  Eqs.|)l5) 
and  (91)  and  the  gradient  operation  of  Eq.  {129): 

SL—OiC.  (r.-S  -  ^ 


The  orthogonality  relation  is  proved  from  Eqs, 
(142)  and  (123) i 


Here  a,  B,  C  are  shown  in  Figure  3  and  is 

given  by 

(144) 

5.6  Eigenfunction  Expansion  of  the  Tensor 
Green's  Function  for  Infinite 
Surface  Recombination 
Velocity 


Having  found  a  complete  set  of  orthogonal 
vector  eigenfunctions  for  the  infinite  surface 


39 


NAVORD  Report  576? 


recombination  velocity  case,  we  construct  a 
tensor  Green's  function.  Affixing  the  super¬ 
script  00  onto  the  parameters  in  Eq.  (126), 

Green's  tensor  is 

rcnn)  .  U^) 

KKlmn  ls  derived  from  Eqs.  (144),  (88),  and  (27): 
\C°  *-(ot)  +(TTI/Q.)JL4(rrm/B)1‘  +(irn/tf+t‘^/D-  (146) 


The  tensor  product  of  a  vector  eigenfunction  pair 
is 


.  To® 

j  ir)  J(ri)  a 


j:Jr« 

jz& 

irk  ‘ 

(147) 

J  yJt. 

jy  Sya 

jyJ*. 

* 

£ 

i'oo  loO 

Jg  jyo 

fee  |00 

JfcJ*  0 

where  the  subscripts  in  the  vector  eigenfunc¬ 
tions  are  not  shown.  The  denominator  of  Eq. 
(145)  can  be  written  as 


H  * 


where  TV  lmn.KKimp,  K  ,  Klm_,  and 
by  Eqs.  (14?) ,  (146) ,  (27),  (144), 
respectively. 

Using  the  abbreviations 


A" 

and 


are  given 
(B?4) 
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cos  ('Prru/e)  *  cr  ;  (rrKu/e)s  sr  ;  (149) 

COS  (.TTru, /&  )  »  Crel  sin(iriru.o/&)=LSr0  , 


where  the  subscript  r  standB  for  1,  m,  or  n 
while  u  stands  for  x,  y,  or  z  respectively. 
Green's  tensor  function  for  the  infinite  sur¬ 
face  recombination  velocity  case  is  given  as 
Equation  (150)  on  page  42. 

When  the  x  dimension  of  the  rectangular 
paralxepiped  goes  to  infinity,  a  technique  sim¬ 
ilar  to  that  used  in  Section  4.3  transforms 
P  00  (r/r0)  into  Eq.  (151)  which  is  given  on 
page  43. 

In  Eq.  (151)  there  are  two  classes  of 
components:  one  contains  ,  obtained  from 
Eq.  (144)  with  t  equal  to  zero,  while  the  other 
contains  KjL,  Eq.  (117).  To  understand  the 
meaning  of  the  two  classes  of  terms,  let  us 
review  the  physical  significance  of  the  Green's 
tensor  function.  Each  tensor  is  composed  of 
nine  components.  When  the  three  diagonal  terms 
are  excited  with  Dirac  delta  functions,  the 
Green's  tensor  describes  the  state  of  the  sys¬ 
tem.  This  state  depends  on  the  properties  and 
the  geometry  of  the  medium  and  the  sources. 

For  the  semiconductor  the  properties  appear  as 
the  diffusion  constant  D  and  the  time  constant 
t.  One  of  the  properties  of  the  source  is  its 
frequency. 

The  components  containing  depend  on 

the  properties  of  the  medium  and  the  frequency  cu 
of  the  sources  qb  well  as  on  the  geometry. 
However,  the  terms  containing  KJJL  depend  only 
on  the  geometry  of  tne  medium  ana  therefore  can 
contribute  only  an  additive  constant  to  our 
final  result.  We  have  examined  thi3  constant 
and  found  it  without  physical  meaning. 


41 


NAVORD  Report  5762 


c_ 

42 


Equation  (150) 


// 
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Iquation  (151) 

[ote:  When  two  signs  are  given  the  upper  one  is  used  when  x>x 
and  the  lower  one  when  x  <x. 
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Therefore  to  save  space  we  neglect  the  terms 
depending  on  in  our  further  analysis; 

In  the  sud sequent  work  xfl>x  and  the  upper 
signs  are  used  in  Sq.  (151).  The  tensor  Green's 
function  for  s  - and  a=“°  is  given  by  Eq,  (152) 
on  page  45.  Equation  (152)  is  the  required 
expression  for  the  tensor  Green's  function  for 
the  three  dimensional  p-n  junction  with  infinite 
s  on  the  transverse  faces. 
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CHAPTER  VI 

NOISE  CURRENT  SPECTRUM  IN  THE  p-M  JUNCTION 
WITH  INFINITE  SURFACE  RECOMBINATION 
VELOCITY 


6.1  Introduction 

We  have  derived  explicit  expressions  for 
the  noise  generators  appearing  in  the  Langevin 
equations  of  the  semiconductor,  and  for  the 
scalar  and  tensor  Green's  functions.  We  now 
use  Green's  functions  to  stun  the  contributions 
from  the  infinitesimal  noise  sources  and  derive 
expressions  for  the  total  noise  spectrum  of  the 
p-n  Junction.  The  recombination  noise  spectrum 
is  derived  from  the  scalar  hole  density,  Eq. 

(82),  and  the  scalar  Green's  function,  Eq.  (114). 
The  diffusion  noise  spectrum  Is  found  from  the 
vector  diffusion  current  density,  Eq.  (121), 
and  the  tensor  Green's  function,  Eq.  (152). 

Since  the  recombination  and  diffusion 
noise  spectra  are  the  result  of  independent 
elementary  processes,  they  are  derived  separ¬ 
ately.  These  noise  spectra  are  added  together 
to  obtain  the  total  noise  spectrum. 

We  consider  first  the  case  of  infinite 
surface  recombination  velocity  because  we  have 
been  able  to  derive  both  the  scalar  and  tensor 
Green's  functions  for  this  case. 


6.2  Recombination  Noise  Current  Speotrum 

Let  us  examine  the  surface  integrals  in 
the  expression  for  the  hole  density,  Eq.  (82): 
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p *°(r)ls  -flG“*V0  ptrj- J5lrc)  Ve  Ga°]-dAe  .  (153) 

The  subscript  s  denotes  the  contribution  from 
the  surface  integrals,  and  the  superscript  oo 
indicates  that  we  are  considering  the  case  of 
infinite  surface  recombination  velocity.  Since 
the  excess  hole  density  p°°(r0)  on  all  the  sur¬ 
faces  equals  zero  for  the  boundary  conditions 
discussed  in  Section  4.4,  the  second  term  on  the 
right  hand  side  of  Eq.  (153)  Is  zero.  Further¬ 
more  an  examination  of  Eqs.  (114)  and  (115) 
shows  that  Q°°  vanishes  at  the  boundaries  Blnce 


and  F.  ,<r0) 

n)  n 


20 


Hi*8 

b*  *  ® 


|x*x4|20(154) 

CM.  =A"-'(155) 

'to-O 


In  each  of  these  expressions  there  are  two 
evaluation  surfaces;  each  expression  is  to  be 
•evaluated  separately  at  each  surface.  There¬ 
fore,  the  first  term  on  the  right  hand  side  of 
Eq.  (153)  is  zero  and  there  is  no  contribution 
from  the  surfaces.  Thus  the  excess  hole  den¬ 
sity  equation  (82)  becomes 

| c~(r)=J(*  Sdva.  (156) 

The  excess  hole  density  is  transformed  to 
current  density  at  the  x*0  plane  with  the 
homogeneous  form  of  Eq.  (14;,  and  the  current 
density  is  integrated  over  the  xaO  plane.  The 
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JtM  I  fS(rt)dtJ%<h0j  J *G?dx  Dv  (157) 

Since  S(rQ)  is  not  a  function  of  the  observa¬ 
tion  coordinates,  wo  get 

f (aa/^i  ’(i58) 


where  Q 00  is  given  by  Eqs.  (114),  (115) >  and 
(116).  Substituting  Eq.  (158)  into  (157)  the 
current  at  the  xaO  plane  is 


Ir(f>  =  U59) 


-d,  («./r‘)2  rnxic  >  «?<-*•  . 

55  ©  0 

From  now  on  the  fact  that  the  current  is  in  the 
x-direction  and  is  evaluated  at  the  x«0  plane 
will  not  be  shown.  The  subscript  r  designates 
the  recombination  process. 

A  noise  process  is  characterized  by  its 
spectrum.  To  derive  the  recombination  current 
spectrum,  w(U,.  <♦)!*’  ),  we  use  the  relation 
between  the  spectrum  and .the  Fourier  frequency 
components  given  by  Rlge1 ' : 

u/(ji/)  =  L«m  l,(  *)!*(  t)/T,  (160) 

r*°° 

where  T  is  the  time  interval.  Substituting  Eq. 
(159)  into  (160),  using  the  relation  for  the 
recombination  source  spectrum. 
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T-*  ’** 

and  employing  ff<t)6(X-*o)dx  (Vo)  t  (162) 


the  recombination  spectrum  becomes 


ur(|rr|i)  = 


(163) 


rrf  a' 

®*4  p 

The  symbols  w(|Srr),  smA,  and  Kumn,  are  given 
by  Eqs.  (6  5)  ^  (l?9)j  anc?°(ll7).  fRe  Dirac  delta 
function  in  Eq.  (161)  expresses  the  fact  that 
at  each  point  of  the  semiconductor  it  is  as¬ 
sumed  that  the  noise  source  is  uncorrelated 
with  the  noise  sources  at  all  other  points. 

The  recombination  source  function,  Eq. 
(65),  contains  <J?t(r0)>  ■  Thia  quantity  is 
defined  by  Eqs.  (fill)  to  (B13)  and  is  made  up 
of  two  parts,  the  thermal  equilibrium  hole  den¬ 
sity  <pfc(r0)/Nt£md  the  average  excess  hole 
density  <pt(r0)>E,  Integrating  the  part  of  Eq. 
(163)  which  contains  the  thermal  equilibrium 
hole  density,  the  Nyquist  current  spectrum  is 
obtained : 


ur(lbl/)„=fl'6  VOCl]  *.  (164) 


where  A°°  and  Ki/mn  are  given  by  Eqs.  (B34)  and 
(117)  respectively. 

The  excess  recombination  noise  spectrum 
is  derived  by  integrating  the  part  of  Eq.  (lb}) 
which  contains  the  average  excess  hole  density, 
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Eq.  (B21).  Before  integration  thia  equation  is 


u No) 


(165) 


y  f  f  [  I  '7 , ‘a  ^  •><>  +  ^ 0  ' 

o  J«  to  to' tv"  m  n'n" 


Integrating  Eq.  (165)  yields 

urOtfOg  *[$#]  ¥ 


(166) 


H  n'n" 

odd 

where  the 
ard  given 


\mn  +  ■* 


symbols  A 80  .  Pb(0),Smni  *Cmn,  and  KU»mn 
by  Eqs.  (B?4),  < ^$)?°(l4§^  117)  tnST 


ard  given  by  Eqs.  (B54),(B$),  (l4§7,  (117) 
(B22;  respectively.  In  deriving  Eq.  (166) 
following  Integral  has  been  used: 


(*r A'o  W «.  3  -4fl r  r'r"l4 r*  *  *r‘ Y- tl  r'r*t\\  ( 167 ) 


where  only  those  values  of  r,  r',  r"  appear  for 
which  the  value  of  |r*r*r"|  is  odd.  The  symbol 
r  stands  for  m  or  nj  9  for  B  or  C;  and  u  for  y 
or  z,  respectively.  Therefore  in  Eq.  (166) 
there  appear  only  those  values  of  m,  m',  m"  and 
n,  n‘,  n'1  for  which  the  expressions  |m*m±m“|  and 
|n*nW|  are  odd. 
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6.3  Diffusion  Noise  Current  Spectrum  for 
Infinite  Surface  Recombination  Velocity 

The  expression  for  the  diffusion  current 
density  is  given^by  Eq.  (121)  with  <»  super¬ 
scripts.  Since  $']“and  v.-f^go  to  zero  at  the 
boundary,  there  is  no  contribution  from  the 
boundaries.  For  the  x  component  of  the  current 
density,  Eq.  (121)  becomes 


jsw * f j'C* + nir,*  +•  a v0 .  (i68) 


To  get  the  total  diffusion  current  IQ(f)  in  the 
x  direction  at  the  x=0  plane,  Eq.  (lo8)  is 
integrated  over  this  plane.  The  result  is 


dj<M*i(l69) 


Performing  the  integration  az  the  x»0  plane, 
using  Eq.  (152),  the  diffusion  current  becomes 


(170) 


i.U>' ft 


Tf 


5. 


K 


w,Vi 


■m  y\ 


v\B  tn  C  / 


where  Smo 
are  the  u 
The 


and  Cmo  are  given  by  Eq.  (149)  and  SDu 
components  of  the  diffusion  sources, 
diffusion  current  spectrum  is 
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W(!I^-  (171) 


Here  we  have  used  an  equation  similar  to  (160) 
and  the  diffusion  source  spectrum  given  by 


Hi  S.?>^  a  wQ5b  J]6U-r.)  dv.  .  (172) 

At  an  arbitrary  point  the  diffusion  sources  in 
the  three  orthogonal  directions  are  independent. 
Furthermore,  the  diffusion  sources  at  an  arbit¬ 
rary  point  are  assumed  to  be  independent  of  all 
other  points. 

In  Eq.  (171)  the  noise  sources,  Eq,  (75)* 
contain  the  average  hole  density  <p£V)>.  As  in 
Section  6.2  the  part  of  Eq.  (171)  which  contains 
the  thermal  equilibrium  hole  density  gives  the 
Nyqulst  noise  while  the  part  which  contains  the 
excess  hole  density  gives  the  excess  noise. 

Integrating  the  Nyqulst  part  of  Eq.  (171), 
the  Nyqulst  diffusion  noise  is 


fitter* 


(173) 


where  A°°  and  are  from  Eqs.  (B^4)  and  (117) 
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respectively,  and  is  from  Eq.  (144)  with  £ 
equal  to  zero. 

The  excess  noise  part  of  Eq.  (171)  is 


w<WV  (A.V’,‘Dpk(<i/ir6) 


{•■pUCC.. + k;;v  +  K- dv.} 


■\T\  YaSvi*'  Vl  w 1  y\" 

I  ^tw&Cwa  5tm*ri  5y\<i  Syi‘o  5y\‘'c 

-rr»"  vvw 

^  Sw*  t-r.C-a'j,  5„-|, 

KVTI^t*'  t\’‘  C  *■ 


(174) 


vn  tn’tw*' 
y\ 

o<i«( 


Here  Eq.  (B21)  is  substituted  for<pt<0>in  (75), 
which  in  turn  is  put  into  (171). 

Integrating  Eq.  (174)  shows  the  excess  dif¬ 
fusion  current  spectrum  to  be 


(175) 


Zj,  K' "  w« '  +  ^v?«V'  IP-"1'’  W 1  -  -m"  *•)  -  (t»V)  w  *-  tH 


,vn 

oA<\ 


where  A  ,  PbPjKKmn*  and  K^inm  are  given  by  Eqs. 
(BJ4),  (B5),  (117),  and  (B22) .  Here  we  used  Eq. 
(167)  and  the  Integral 
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[ - (^r 0/rt)(r^ r-*- v-^V'-r^Vryff,  (176) 


where  r  is  either  m  or  n.  In  Eqs.  (167)  and 
(176)  only  those  values  of  r.  r*,  rft  appear  for 
which  the  expression  |rtr,tr',|ls  odd. 


6.4  Total  Nyquist  Noise:  Stochastic  Theory 
and  Nyquist  Law 

The  total  noise  spectrum  calculated  from 
stochastic  theory  is  composed  of  two  parts: 
the  recombination  current  spectrum  and  the  dif¬ 
fusion  current  spectrum.  Since  the  recombina¬ 
tion  and  diffusion  processes  are  independent 
statistically,  we  add  the  spectral  densities. 
For  the  total  Nyquist  noise  current  spectrum 
Eqs.  (l64)  and  (173)  are  summed  to  get 

[>[ 

H1T| 

The  symbol  ft.  (KuLd)  means  the  real  part  of 
From  equations  |e$0)>  $322)  and  (117;  we  get 

fi(K^=k  te. *;Ut  ■  n*3) 


Z  (L(,¥l 


o^rrj 


(177) 


With  Eq.  (178),  Eq.  (177)  becomes 


04* 


(179) 


As  a  check  on  the  above  method,  the  Nyquist 
noise  current  spectrum  is  calculated  from  the 
conductance  at  the  x=0  plane.  From  the  Nyquist 
law1^  the  noise  spectrum  is 
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wl\Xr\)M-4-KTG  ,  (180) 


where  G  is  the  real  part  of  the  thermal  equilib¬ 
rium  admittance  of  the  p-n  junction.  At  thermal 
equilibrium  V0  is  zero  in  the  expression  for  the 
admittance,  Eq.  (B33),  and  Eq.  (lOo)  is  identi¬ 
cal  with  Eq.  (l79).  Thus  the  result  of  the 
stochastic  analysis  is  correct,  giving  us  con¬ 
fidence  in  the  method  of  tensor  Green* a  func¬ 
tions.  We  shall  see  in  Chapter  VII,  on  the 
other  hand,  that  the  stochastic  result  obtained 
by  the  use  of  scalar  Green's  functions  does  not 
check  the  result  obtained  from  the  Nyquist  law. 


6.5  Total  Excess  Current  Spectrum  with 
Infinite  Surface  Recombination 
Velocity 

Adding  the  excess  recombination  and  the 
excess  diffusion  current  spectra,  Eqs.  (165)  and 
(175),  the  total  noise  spectrum  is  found  to  be 


\N 


iu%  -  [i'vbm » 

X  [pt*  t  (rf*  -Ml 


(181) 


are  given  by  Eqs. 

r  represents 


where  A  ,  pp, %mn>  and  Kk*mn  are  B1 
(BJ4),  (B5),  (1177,  and  (B2^.  Here  _ 
either  m  or  n  and  only  those  values  of  r,  *  , 
appear  for  which  |rsr'*r"|is  odd. 

It  is  interesting  to  see  if  the  excess 
spectrum  is  proportional  to  the  steady  current 
flowing  at  xcO.  The  dc  current  density  is 
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derived  from  Eqs.  (12)  and  (B2l).  Integrating 
over  the  x=0  plane  gives  for  total  current 


I“*  P„ID) H°/xC]  5]  /W-W-.  (182) 

Comparing  the  dc  current  with  the  excess 
current  spectrum,  Eq,  (l8l),  the  two  relations 
can  be  related.  However,  the  factor  of  pro¬ 
portionality  is  a  complicated  function  of 
frequency,  semiconductor  parameters,  and  geometry. 


6.6  Divergence  of  the  Expression 
for  the  Nyquist  Noise 

Let  us  examine  the  noise  spectra  for  con¬ 
vergence  at  each  frequency.  The  Nyquist  current 
spectrum,  Eq.  (17?),  is  converted  into  a  double 
integral  and  integrated.  All  terms  in  Eq.  (177) 
are  positive,  and  we  investigate  the  range  for 
large  m  and  n.  Making  the  transformations 


-rrSrv 

B 


•=  x 


A  _  It  Tfri  4  it 

1  =  T  jr 


(18?) 


and  neglecting  1/Dr  and  <o/D  with  respect  to  m 
and  n,  the  series  in  Eq.  (177)  becomes  when  Ax 
and  Ay  go  to  zero 

R  (  f  d*  dj 

+*  '  BC  / 

W\v>l  HV>» 


I 


(184) 


With  the  transformations 


X  =•  j)  cos  0  siflQ  ,  d*4j -j3dj>d0*  (1^5) 
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Eq.  (184)  becomes 


r  i't"1  (\  [ 

I  .  >v\’-  -vv1  Q  C.  I  a  *■  I 

-m.-r,  0.  Jrt 


B.  -H. 
?.  c.. 


4  3 


(186) 


Therefore  w(  |l2|  )N  becomes  infinite. 

To  determine  the  physical  significance  of 
this  divergence  we  investigate  the  admittance 
with  Vp*0  at  the  xaO  plane,  since  the  Nyquist 
current  spectrum  is  proportional  to  its  real 
part.  For  arbitrary  surface  recombination  vel¬ 
ocity  the  admittance  at  the  xsO  plane  is  given 
by  Eq.  (B28).  The  factors  in  the  denominator 
of  the  infinite  series  are  each  bounded  and  are 
neglected  for  discussions  of  convergence,  Fur¬ 
thermore,  from  the  boundary  conditions  Eqs,  (96) 
and  (97) 

,  (187) 

where  r  stands  for  m,  n;  <*  for  b,  c  respectively. 
Only  for  infinite  surface  recombination  velocity 
does  A9r  always  remain  tf/2  independent  of  r. 
Otherwise  A9r  approaches  zero  as  r  goes  to 
infinity. 

If  we  take  r  finite  but  so  large  that 
A9r  is  small,  *>  C(D*)"'  ♦iw/DJ  ,  and 

the  Identity  for  large  positive  integers  m  and 
n 

+  ,  (188) 

the  series  in  Eq,  (B28)  becomes 

Y  ft™  bt  3^  .  (189) 
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As  m  and  n  individually  go  to  infinity  the 
terms  go  to  zero  because  am)*or  (den)*1  goes 
to  zero.  The  value  of  the  input  admittance  is 
finite  except  when  A9p  is  ir/2,  when  it  is 
infinite.  When  saoc  the  surface  of  the  semi¬ 
conductor  seems  to  be  covered  with  a  perfectly 
conducting  layer  which  short-circuits  the  semi¬ 
conductor.  Therefore  the  Nyquist  noise  current 
spectrum  becomes  infinite  when  the  surface 
recombination  velocity  is  infinite. 

A  convergent  expression  for  the  case  of 
arbitrary  s  can  be  written  down  directly  from 
the  Nyquist  Law,  Eq.  (180)  and  the  real  part  of 
the  admittance,  Eq.  (B3I).  This  is  discussed 
further  in  Chapter  VII  for  the  case  of  arbitrary 
s.  Now  we  examine  this  expression  for  the  case 
of  large  but  finite  3. 

From  the  boundary  conditions,  Eqs,  (96) 
and  (97), 


Y  -  s<*/P  =(^a  tan  ft,*. 

When  V  is  large  we  can  write  /9*.  as  ^ 


nr  rir 
z  zy 


rj r 

z. 


tY 


» 


(190) 

(191) 


where  r  is  odd.  Substituting  Eq.  (B31)  into 
(180),  letting  f  be  large  and  the  do  voltage  at 
x»0  be  zero,  we  get  for  the  range  of  small  Ey  , 


tLOSt*)irD/as&)  cos  (nrrZUsc) 
(mir/tX)  -  D/s  &X  n  fr/z)(i  -£>/■*  e;J 


sin  (mrrJD/s  6)1 

hitU+d/s  fi)J 


"  ['  + 


-t 


(192) 


X  |  (l(DrT+[(rnrr//J)(l-'D/st)f'+Ufiir/zc)(i~D/sc)f]l\  tf/tf 

+  (Dt)'\[(.mir/te)(i~j)/j6) f  +[  (mr/zcXi-Mcrffyy 
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If  we  let  sb/D  and  sc/D  become  very  large 
while  £m  and  <f„  remain  small,  Eq.  (192)  becomes 
identical  with  the  Nyquist  current  spectrum  with 
s  s«p  ,  Eq.  (179);  to  make  the  two  equations  agree 
the  conversion  factors  between  Figures  2  and  J> 
are  used: 

B-lb)  02c.  (193) 


6.7  Convergence  of  the  Series  for  the 
Excess  Noise  Spectrum 

To  discuss  the  convergence  of  the  excess 
noise  spectrum,  we  investigate  Eq.  (l8l).  In 
the  denominator  of  this  equation  there  are  two 
factors  which  contain  minus  signs,  one  in  m  and 
one  in  n.  Taking  the  factor  in  m  (identical 
results  are  obtained  with  n).  the  indices  m,  m' 
and  m"  are  related  so  that  must  be  odd. 

For  the  whole  factor  to  be  zero 

( ml  -  m"1)2  (2.  mV  fl  (194) 

Solving  we  get 

7n  ±m' ±m" -  o .  (195) 

Since  zero  is  an  even  number,  the  factors  cannot 
vanish.  Furthermore  the  factor  in  m  is  nega¬ 
tive  whenever  m<mf4m". 

The  question  of  convergence  of  the  excess 
noise  spectrum  is  a  very  difficult  one  to 
answer.  Converting  the  sum  to  an  integral  is 
not  permissible  since  some  of  the  values  of 
r,  r'  and  r"  are  not  present.  Furthermore  only 
conditional  convergence  and  not  absolute  conver¬ 
gence  is  required.  Putting  test  values  into  the 
series  the  denominator  increases  very  much 
faster  than  the  numerator.  It  appears  that  the 
sum  converges  rapidly,  but  a  rigorous  proof 
would  require  numerical  evaluation  of  the  excess 
noise  spectrum. 
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Equation  (l8l)  can  be  used  to  calculate 
the  excess  current  spectrum  for  infinite  sur¬ 
face  recombination  velocity.  Prom  the  above 
discussion  of  the  Nyquist  noise,  we  can  expect 
that  Eq.  (181)  provider  an  upper  bound  for  the 
excess  noise  spectrum  when  the  surface  recom¬ 
bination  velocity  is  large  but  not  infinite. 


6.8  Contribution  of  Electron  Density  Fluctuation 
in  the  p-Type  Material 

Our  analysis  has  only  considered  hole 
conduction  in  the  n-region.  Shockley7  has 
shown  that  electron  conduction  in  the  p-reglon 
is  simply  an  additive  effect  to  the  hole  cur¬ 
rent.  The  noise  resulting  from  the  concentra¬ 
tion  fluctuations  of  electrons  is  statistically 
independent  of  the  hole  fluctuations  and  the 
analysis  is  similar  to  that  made  above.  When 
the  p  and  n  symbols  in  the  expressions  for  the 
hole  fluctuations  are  interchanged  and  when  the 
values  of  D  and  t  for  electrons  in  p-type 
material  are  used,  the  derived  spectra  pertain 
to  electron  fluctuations  in  the  p-type  material. 

The  total  current  spectrum  is  the  sum  of 
the  hole  current  spectrum  and  the  electron  cur¬ 
rent  spectrum. 


6.9  Summary  of  Chapter  VI 

Using  tensor  and  scalar  Green's  functions, 
we  have  derived  explicit  expressions  for  the 
thermal  equilibrium  noise  spectrum,  Eq.  (179), 
and  the  excess  noise  spectrum,  Eq.  (lol),  of  the 
p-n  Junction  when  the  transverse  surface  recom¬ 
bination  velocity  approaches  infinity.  In 
actual  practice  the  surface  recombination 
velocity  is  never  Infinite,  but  may  be  very 
large.  Since  Eq.  (180)  is  convergent,  it 
should  be  used  to  calculate  the  thermal  noise 
when  the  surface  recombination  velocity  is  large, 
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tout  finite.  The  excess  noise  spectrum  can  toe 
evaluated  from  Eq.  (l8l)  and  should  give  a 
good  upper  bound  for  the  noise  when  the  sur¬ 
face  recombination  velocity  is  large,  but 
finite . 
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CHAPTER  VII 

NOISE  IN  A  p-n  JUNCTION  WITH  ARBITRARY 
SURFACE  RECOMBINATION  VELOCITY 


7.1  Introduction 

In  the  discussion  of  the  noise  current 
spectrum  with  an  infinite  surface  recombination 
velocity,  we  had  a  very  important  check  on  our 
work.  We  derived  the  Nyquist  noise  current 
spectrum  by  two  independent  methods:  the  Nyquist 
law  using  the  input  conductance  and  the  stochas¬ 
tic  process  theory  using  basic  physical  prin¬ 
ciples.  These  spectra  are  identical.  In  thiB 
chapter  we  get  the  current  spectrum  for  the  case  'of 
finite  surface  recombination  velocity.  When  we 
compare  the  current  spectrum  of  the  Nyquist 
noise  derived  from  the  Nyquist  law  and  from 
stochastic  theory,  we  find  that  only  for  small 
values  of  surface  recombination  velocity  do 
the  two  methods  agree.  Since  the  stochastic 
process  method  uses  basic  physical  principles 
and  rigorous  methods,  we  present  it,  but  we  do 
not  resolve  the  question  of  which  result  is 
correct. 

For  the  case  of  arbitrary  surface  recom¬ 
bination  velocity  we  use  the  rectangular  coor¬ 
dinate  system  shown  in  Figure  2  with  the  p-n 
junction  at  the  x*0  plane.  The  origin  is  at 
the  center  of  the  p-n  junction. 

To  find  the  current  spectrum  for  the 
three  dimensional  p-n  Junction,  the  excess 
hole  density  p  is  derived  with  the  aid  of  the 
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scalar  Green1 3  function.  This  density  p  is 
obtained  from  Eq.  (82).  Since  the  Green's  func¬ 
tion  satisfies  the  same  homogeneous  boundary 
conditions  as  the  unknown  function,  Eqs.  (96) 
and  (97) ,  the  surface  integrals  in  Eq.  (82)  are 
zero.  For  the  surfaces  xo=0  and  x0= 00  this  is 
proved  with  a  technique  similar  to  that  used  in 
Section  6.2. 

The  vanishing  of  the  remaining  surface 
integrals  is  proved  in  the  following  manner.' 

They  can  be  written  as 

pwl,  =/IG  ^  -  rCr*>  fl. .  .  (196) 


and  the  mixed  boundary  conditions,  Eq.  (93), 
can  be  written  as 


»*(fy+i)«p/D,  (197) 

where  the  symbol  ^  is  defined  by 


JL 


1 


■» 


(198) 


Green's  function  satisfies  Eq.  (197).  To 
prove  this,  Green's  function,  Eq.  (109 ) ,  is 
written  as 


(199) 


and  is  differentiated  to  get 
jG  _ r-AM1*.  ] o 


(200) 


where  r  denotes  m  or  n,  f  denotes  ju.  or  u  ,  and 
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u0  denotes  y0  or  zQ  respectively.  Eqs,  (199) 
and  (200)  are  evaluated  at  the  boundaries  with 
the  aid  of  the  eigenfunction  form  of  the  mixed 
boundary  conditions,  Eqs.  (96)  and  (97).  When 
we  compare  the  resulting  equations,  we  see  that 
Green's  function  satisfies  Eq,  (197). 

Substituting  Eq.  (197)  f-r  N£p  and  a  simi¬ 
lar  expression  for  G  into  Eq  (196),  p(h)|s  is 
found  to  be  identically  zero  and  Eq.  (82)  becomes 

p(r)=yS(OGO,lO  <Av0  .  (201) 


Here  Scfcjand  G  are  given  by  Eqs.  (35)  and  (109) 
respectively. 


7.2  Bulk  Recombination  Current  Spectrum 
for  Finite  Surface  Recombination 
Velocity 

The  recombination  current  spectrum  deriv¬ 
ation  is  exactly  the  same  as  that  for  the 
infinite  surface  recombination  velocity.  The 
Nyqulst  current  spectrum  is 


»0U‘)  ,  <S0£?> 


where  A,  Lmn,  and  Kymn  are  given  by  Eqs. 

(B20),  (B19T  (lSSJ  and  (ffg) j  and  the  symbols 
q  and  t  are  defined  after  equation  (12). 

The  excess  current  spectrum  is 

^(IU’)c  = 

V’ <”•  «'"■  *-V  *  «, Jf 05 ’ 

n,  m*(  nM 


64 


NAVORD  Report  57^2 


where 
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(204) 


E“  'L c"  c’’- du- '  i<x*  «av;J 

X  i^N  ’■A-  «i) 
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f 

cro* 

sro>  cr<x  > 
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*  W>jS,  , 

s.t  *  u. , 

c*< 

5  , 

Sr«  *  *,h/?ra<  ' 

(205) 


where  r  stands  for  m  or  n;  at  for  b  or  c;  and  uQ 
for  yQ  or  zQ  respectively.  pi,(0)  13  given  by  Eq. 
(B5) .  The  symbol  Ki£imn  Is  given  by  Eq.  (B18) 
and  the  other  symbols  are  Identified  after  Eq. 
(202 ). 


7.3  The  Diffusion  Current  Spectrum  for  Finite 
Surface  Recombination  Velocity 

The  diffusion  current  spectrum  for  the 
finite  transverse  surface  recombination  velocity 
requires  a  new  approach.  We  reexamine  the  solu¬ 
tion  for  the  excess  hole  density,  Eq.  (201). 

By  transforming  this  equation  the  diffusion 
noise  source  takes  on  the  same  form  as  that 
derived  in  Chapter  III.  Gradients  of  the 
scalar  Green* s  function  result  from  this  trans¬ 
formation. 

We  start  the  derivation  of  the  diffusion 
current  spectrum  with  the  solution  for  excess 
hole  density.  The  source  function  S  is  now  Sp, 
Eq.  (35b  and  Eq.  (201)  becomes 

PD<r>=/[V.V0;>]G«Wo 
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The  appearance  of  the  Lapiacian  in  the  source 
function  makes  it  impossible  to  use  the  physical 
source  function,  Eq.  (75).  We  transform  Eq. 
(206)  into 


p(r)  *  j [Q  r.  Po]  -art.-/ [v.  Pi]  •  [v.  G]  av. 


(207) 


with  the  vector  identity 

G  [v/p]  -  v«  '[Gtip]  -  tep]*  [V.o]  (208) 

and  Gauss*  theorem. 

The  surface  Integrals  in  Eq.  (207)  result 
from  the  diffusion  current  to  the  transverse 
surfaces  of  the  semiconductor.  This  current  is 
caused  by  the  surface  recombination  and  is 
related  to  the  surface  recombination  velocity. 

It  is  important  to  note  that  the  surface  sources 
did  not  enter  prior  to  the  transformation  of  Eq. 
(207).  The  three  volume  integrals  are  due  to 
the  volume  diffusion. 

In  order  to  derive  the  diffusion  current 
spectrum,  we  solve  for  the  total  diffusion  cur¬ 
rent,  ID(f),  in  the  x  direction  at  the  x»0  plane. 
When  Eq.  (20/)  is  written  in  component  form  and 
integrated  over  the  x=0  plane,  the  diffusion 
current  becomes 

ID«  3  ^  )  ***•  (  209  ) 

x  (-/=i 1.)  j  «.  <i?u' j  -fti  c  (•ifi') 

■{£  4  c  &|  ja?. 
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Here  pQ  is  substituted  for  p(r0).  The  symbols 
A,  KfeJ  and  p  are  glv§n  by  Eqs.  (B20), 

(B19),  (II277  and  (96;  respectively.  The  sym¬ 
bols  D  and  q  are  defined  after  Eg,  (12).  In 
deriving  Eq.  (209)  the  following  facts  and  equa¬ 
tions  are  used:  the  scalar  Green's  function  is 
zero  at  the  planes  x=>0  and  x  =  «»>  Eq.  (14) j  in 
Eq.  (207)  only  Green's  function  depends  on  the 
observation  coordinates. 

The  diffusion  current  spectrum  is  derived 
by  substituting  the  diffusion  current,  Eq.  (209), 
into  an  equation  similar  to  (160)  and  is 


W)  -4D  (210) 

x  [o1!/  <pt<b)>  cJb  i  «&  -^<pt(c)>  c~* 

+l  L <f5t>  c"#  c^c-°  c*‘  ^ 

+aa  C  sl  ]  «**/•**  *•]  J  * 


where  the  symbols  cL  sJL, 
defined  by  Eq.  (205;,  while 


,  and  s^.^  are 
the  other  symbols 

are  Identified  after  Eq.  (209).  The  factor  two 
in  the  surface  integrals  comes  from  the  alge¬ 
braic  addition  of  the  two  equal  uncorrelated 
noise  spectra  produoed  by  currents  to  opposite 
transverse  surfaces.  In  the  derivation  of  Eq. 
(210)  the  relation  characterizing  the  surface 
sources  is 


TLr -f  ^ l  ■  <211> 
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while  the  relation  characterizing  the  bulk 
diffusion  sources  is 


Jr'"1  T[(t£Xi£f]  •  (212) 

In  Eg.  (211)  uQ  stands  for  y0  or  z0  while  in 
Eq.  (212)  u0  stands  for  x0,  y0,  or  zQ;  /„  denotes 
any  coordinate  except  u  .  In  Eq.  (211)  the 
symbol  ws  ( |f£|J  is  gi9en  by  Eq.  (80)  while  in 
Eq.  (2127  the  symbol  wD(||fcr)  ls  derived  from 
Eqs.  (7?),  (6o),  and  the  homogeneous  form  of  (14). 

To  evaluate  the  diffusion  current  spectrum 
it  is  divided  into  Nyqul3t  noise  and  excess 
noise.  For  the  Nyquist  noise  p  is  substituted 
for  ^Pt(r)/>  in  Eq.  (210)  and  the  equation 
integrated,  (213 ) 

win/)  *  Vfl*Pw!iYV!!&- 

+I  k^c!: 

m.n.W  mln 


1*3?  * 

where  the  integrals  L  ,  Lr,  and  Lr  are  given 
by  Eqs.  (l?2),(  (103),  and  (133)  respectively; 
when  r=r‘  ,  Lrr'  is  replaced  by  s  ls  the 
surface  recombination  velocity.  The  other 
symbols  are  identified  after  Eq.  (209) . 

For  the  excess  noiBe  spectrum,  p^r),  Eq. 
(B17),  is  substituted  for<pt(r)>in  Eq.  (210). 
Integrating  the  resulting  equation,  the  excess 
diffusion  current  spectrum  is 
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W(II ol%  - 


(214) 


xS  . 

irnv'm  , 
n  n'  n 


1  r-lm)  ~  (n) 

'mn  ^ccc  t-ccc 

^Kmn  +^Kto/n'  +  ^K'tn"n 


r  *  o  Flm)  flfl  EjL 

X  KmnW+^'?*f  +^"'  Ecc'c 


L  Etce  ^tc 


4-  &A 
D  L 

The  symbols  E^„  and  p,  (o)are  given  by  Eqs. 

(204)  and  (B57T^BymbolD  E<:’  is  given  by 

oSC 

r- 

OL 

E«c  “It  5in^o  C05pr„u.  duc=  ^  ^(215) 

=  MPr  ^kl  CKa  v  l‘  C«  6r'«  cr1d  “  Chl  Cri*  SV&) 

~  (£ '£)  (l3,  cw  4^  *  ?•  5r«  CH*  Cr'b»  +#»  5r4  r  <*)] 


f  f 

where  C  ^  and  SrtX  are  defined  by  equation 
(205);  the  other  symbols  are  identified  after 


equation  (209). 


7,4  Total  Nyqulst  Noise  for  Finite  Surface 
Recombination  Velocity:  Stochastic 
Process  Method  and  Nyquist 
Law  Method 

The  total  Nyqulst  current  spectrum  deter¬ 
mined  by  the  stochastic  process  method  equals 
the' algebraic  sum  of  recombination  current 
spectrum,  Eq.  (202),  and  the  diffusion  noise 
spectrum,  Eq.  (213),  and  is 
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m  n 


(216) 


X[Lmn  (  <»*>*'  -Ku/) 

L"i}- 

n*n  The  Nyquist  current  spectrum  for  the  p-n 
Junction  is  now  derived  from  the  Nyquist  law, 

Eq,  (l80),  and  the  input  conductance,  Eq.  (BJl); 

wfl!1!)^  s£Y2TA  -25^^  U 


~mn 


xjti<L*^‘+W* 


(217) 


Comparing  Eqs.  (216)  and  (217)  we  see 
that  the  two  current  spectra  are  not  the  same. 

The  former  has  a  triple  Infinite  series  while 
the  latter  has  a  double  infinite  series.  When 
the  parameters  sb/D  and  sc/D  become  very  small, 
agreement  is  attained  between  the  two  expres¬ 
sions  for  the  Nyquist  current  spectrum.  In  this 
case  the  X  factors  converge  very  rapidly,  as 
shown  by  Shockley, 7  and  the  first  terms  of  the 
series  represent  the  Nyquist  spectrum  adequately. 
For  the  usual  p-n  material  1/Dt  is  large  with 
respect  to  the  p0  values  in  the  Kkoo  factor, 

Eq .  ( 112 ) . 

For  small  scx/d  values  the  current  spect¬ 
rum  by  both  the  stochastic  process  and  the 
Nyquist  law  methods  is 

w  (II'I)'  ■  w(U*l)‘  *4^ D/?[Kk  J  (218) 
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/?[Kj  *2Afl[  [<'/Dr>  *  (219) 

+  (wD"')*]  1  +  £  (DT)"’+  sTf  (*■  *  c 

Concerning  the  difference  between  the 
stochastic  process  spectrum  and  the  Nyqulst 
law  spectrum  for  large  values  of  aa/D,  we  note 
that  the  transformation  giving  Eq.  (207), 
which  was  made  in  order  to  use  the  physical 
diffusion  noise  source,  results  in  the  gradient 
operation  on.  the  scalar  Green's  function. 

While  this  seems  correct  formally,  it  may  be 
that  it  throws  the  problem  into  the  domain  of 
the  tensor  Green* s  function.  The  original 
reason  for  Investigating  the  tensor  Green's 
function  was  to  use  the  diffusion  noise  sources 
without  an  additional  transformation.  Unfor¬ 
tunately  an  explicit  expression  for  the  tensor 
Green's  function  for  the  case  of  arbitrary 
surface  recombination  velocity  has  not  been 
found.  Therefore,  we  are  not  able  to  determine 
whether  the  scalar  method  for  the  diffusion 
spectrum  Is  in  error  or  whether  the  Nyqui3t 
theorem  is  not  general  enough  to  handle  noise 
in  a  p-n  junction  with  mixed  boundary  condi¬ 
tions. 


7.5  Total  Excess  Current  Spectrum  for 
Finite  Surface  Recombination 
Velocity 

The  total  excess  current  spectrum  for 
arbitrary  surface  recombination  velocity  is 
obtained  by  adding  the  excess  current  spectrum 
for  bulk  recombination  and  for  diffusion,  Eqs. 
(203)  and  (214)  respectively.  These  two  sta¬ 
tistically  Independent  spectra  add  algebraically, 
and  the  total  excess  current  spectrum  is 
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(220) 


X 

X 


JJU  mV  tec  Ectc  (K(<rnn 

'prefer 


Hcoi  3  ijMJ. 4au/3, 

Let  us  relate  the  dc  current  in  the  x 
direction  at  the  x=0  plane  to  the  excess  cur¬ 
rent  spectrum.  The  d if fusion  current  density 
in  the  x  direction  is  derived  from  Eq.  (B17) 
and  the  homogeneous  Eq.  (14).  Integrating 
this  current  density  over  the  x=0  plane,  the 
total  dc  current,  Idc,  is 

tc  “Pf^  (2al) 

s6)^(8K)!qB^,ona 

spectively;  the  sym¬ 
bols  D  and  q  are  defined  after  Eq.  (12). 

Comparing  this  dc  current  with  the 
excess  current  spectrum,  Eq.  (220),  we  find  a 
relation  between  the  quantities,  but  the  factor 
of  proportionality  is  a  complicated  function  of 
frequency,  geometry,  and  semiconductor  constants. 
When  s  is  small  the  excess  noise  spectrum,  Eq. 
(220),  becomes 

wUtVg  =  4 gA -  “V&f  J  (sss) 


where  the  symbols  A, 
fan  are  given  by  Eqs.  (B 
(96),  (B5),  and  (111)  re 


72 


NAVORD  Report  5762 

and  the  dc  current,  Eq.  (221),  becomes 

<a»> 

where  Kt.  -  (Dt)  '  •*  i  W'.c'VD.  (S24) 

n  00 

When 


‘'/D  «  > 

W(II‘/)E  »  /i/V  fr  X*  .  (225) 

This  has  the  familiar  form  of  the  shot-effept 
phenomenon  in  a  temperature-limited  diode, ■‘•5 
except  that  the  constant  for  the  diode  is  two, 
while  Eq.  (225)  gives  8/5. 


7.6  Spectrum  of  Electron  Density  Fluctuations 
in  p-Type  Material 

As  in  Section  6,8  the  noise  current 
spectrum  of  the  electron  density  fluctuations 
in  a  p-type  material  is  similar  to  the  spectrum 
of  the  hole  density  fluctuations  in  n-type 
material.  The  derivation  of  the  electron  dens¬ 
ity  fluctuation  spectrum  is  the  same  as  that 
described  in  Section  6.8.  The  total  current 
spectrum  is  again  equal  to  the  sum  of  the  hole 
and  electron  current  spectra. 


7.7  Summary  of  Chapter  VII 

Using  scalar  Green's  functions  we  have 
derived  explicit  expressions  for  the  Nyquist 
noise,  Eq,  (216),  and  the  excess  noise,  Eq.  (220), 
for  the  case  of  arbitrary  surface  recombination. 
Thlq  would  appear  to  constitute  a  complete 
solution  to  the  three-dimensional  p-n  junction 
noise  problem.  However,  because  of  the  lack  of 
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agreement  between  our  result  for  the  thermal 
equilibrium  noise,  Eq.  (216),  and  the  result 
obtained  from  the  Nyquist  theorem,  Eq.  (217), 
it  is  recommended  that  further  study  be  made 
before  complete  reliance  is  placed  on  the 
stochastic  results. 
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CHAPTER  VIII 


CURRENT  SPECTRUM  IN  A  TWO-DIMENSIONAL 
SEMICONDUCTOR  WITH  A  p-n  JUNCTION 


8.1  Introduction 

In  many  noise  experiments  the  semicon¬ 
ductor  sample  has  a  two-dimensional  nature. 

The  sample  is  several  diffusion  lengths  in  the 
longitudinal  or  x  dimension  and  for  all  practi¬ 
cal  purposes  may  be  considered  Infinite  in 
this  direction.  In  one  of  the  transverse 
directions,  say  the  y  direction,  the  filament 
is  narrow  and  the  surface  recombination  velocity 
is  important.  In  the  other  transverse  direction, 
the  filament  is  so  wide  that  the  effect  of  sur¬ 
face  recombination  velocity  is  negligible.  A 
bias  voltage  across  the  p-n  Junction  causes  a 
current  to  flow  in  the  x  direction. 

Since  there  is  no  variation  in  the  z 
direction,  the  del-operator  in  Eqs.  (34)  to 
(37)  becomes 

^  3  (226) 

Except  for  this  change  the  current  spectrum  is 
derived  with  the  technique  used  to  determine 
the  spectrum  of  the  three  dimensional  semicon¬ 
ductor  containing  a  p-n  junction. 
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8.2  Two  Dimensional  Current  {spectrum  for 
Infinite  Surface  Recombination  Velocity 


The  Nyquist  noise  current  spectrum  is 
ur(|Ia)N  =/T(^r)^afVi  ®  £  ft  (***>)/*"' , 


(227) 


where  (228) 

A(\CKr)  «  /l^[([n«/ts]V/ot-)a’+(w^*ji4-([rr«/B]+ 


The  Nyquist  noise  spectrum  derived  with  the 
Nyquist  law  and  the  conductance  of  the  two 
dimensional  p-n  junction  is  the  same  as  Eq. 
(227). 

The  excess  noise  spectrum  is 


£  1  ,  ,  (229) 


where  the  only  values  of  m,  m",  m"  appearing 
are  those  for  which  Imim'-m"!  is  odd.  The 
symbols  pb(0)  and  A“  are  given  by  Eqs.  (B5) 
and  (B34),  respectively;  the  symbols  D,  q  and 
t  are  defined  after  Eq.  (12).  The  symbols 
and  K^im  are  defined  as 


5  [(nm/8)aH(£.u,/p)-M/Dir]4  (2?0) 

K"m  5  (251) 


In  two  and  three  dimensions  the  current 
spectra  behave  the  same  (see  Section  6.6);  the 
Nyquist  current  spectrum  diverges,  but  the 
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excess  current  spectrum  appears  to  converge 
rapidly.  The  excess  current  spectrum  can  be 
used  as  an  upper  bound  for  large  values  of 
surface  recombination  velocity.  As  in  the 
three  dimensional  case,  the  proportionality 
between  the  excess  noise  spectrum,  Eq.  (229), 
and  the  dc  current  at  the  x=0  plane 


T»  _  1%  Dph(o)  fl  “ 
s  ir* 


is  very  complicated. 


(252) 


8.3  Two  Dimensional  Current  Spectrum  with 
Finite  Surface  Recombination  Velocity 

The  Nyquist  current  spectrum  obtained  by 
the  stochastic  process  method  is 


ur(,n)’  *  VflV-0  (255) 

v  m  Kskj 

+  H  JC.K,,(asD'l4 ^/(Kk^K k^1)] 

m  y  *n'  J  * 

The  Nyquist  current  spectrum  obtained  by  the 
Nyquist  law  and  the  input  conductance  method  is 


(234) 


respectively;  and  the  following  expressions 
define  the  remaining  parameters: 


(235) 


given 

(96) 


urt, ZlltbfyZlCLlRlKKj. 

mu. _ _ u. .1 .  .  r  “  r mm ®  _ _ a  Q 


The  symbols  A,  LI, 
by  Eqs.  (B20),  (10 


03)™  ’ 


L" 

(15?), 


(152)? 


are 

and 
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kk»«=  t>n  +  A"  ’  (236) 

*k-*kw«+%-  .  (237) 

+  .  (238) 

When  sb/D  becomes  small,  the  Nyquist 
spectrum  derived  by  the  Nyquist  law  .and  by  the 
stochastic  process  method  is 

♦■*•**, .tJiW  .  <«9) 

where  MO  + js')1*  o']*  .  (£i|0) 


In  two  dimensions  as  in  three  dimensions  the 
Nyquist  current  spectra  obtained  by  the  sto¬ 
chastic  processes  method  and  by  the  Nyquist  law 
method  become  the  same  only  when  sb/D  is  small. 

The  excess  current  spectrum  with  finite 
surface  recombination  velocity  is 

»«tv  Ik. -E™ (V<«J  '2l,1) 

x®>tr+K>.K;.,+A .jwe£(Kf 

+  la  (DEcTc)'  c«»  J5»b  c.t tes^.fc]^  • 


The  symbols  and  Ei£’  are  given  by  Eqs. 

(B5),  (204),  and  respictively;  the  sym¬ 

bols  D,  q,  and  t  are  defined  after  equation  (12) ; 
and  the  rest  of  the  symbols  are  defined  after 
Eq.  (234). 
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The  dc  current  is 


i* V  ^ 

As  In  three  dimensions.  Section  7.5,  a 
complicated  proportionality  factor  relates  the 
excess  current  spectrum  and  the  dc  current  for 
arbitrary  surface  recombination  velocity.  For 
small  sb/D  the  excess  current  spectrum  is 

(243) 

+  (L[(Drr'+5(Dtr,J2'  +  ^I>'')  I  j 

*{  { [[[(PtA 

+  (et>\  suvy'}k + [(Dtr‘+5  uD)'Y/"f  • 

When  sb/D  is  small  and  when  (  'x>/D)«  K^i  ,  the 
excess  noise  spectrum  becomes  '  A  ' 

w(II \  Kk.0.  (244) 

For  small  sb/D  the  dc  current,  Eq.  (242), 
becomes 


Ldc  (2^5) 

Thus  the  relation  between  the  excess  current 
spectrum  and  the  dc  current  is 

V/01%  -(!/*)%  Zdc.  (246) 

The  contribution  from  the  electron  den¬ 
sity  fluctuations  in  p-type  material  is  similar 
to  the  three  dimensional  case,  Section  6.8. 
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8.4  Comparison  with  Experiment 

Experiments  designed  to  study  surface 
phenomena  are  generally  performed  with  thin 
slices  of  germanium  or  silicon  and  fit  the  two 
dimensional  geometry  quite  accurately.  Experi¬ 
mental  methods  are  available  to  vary  the  sur¬ 
face  recombination  velocity.  Certain  gaseous 
ambients,  such  as  water  vapor,  have  been  found, 
to  Increase  the  surface  recombination  velocity^5 
grea.tly.  Our  results  can  be  used  to  predict 
the  effect  of  such  surface  treatment  on  the 
minority  carrier  noise  of  a  p-n  junction.  For 
the  case  of  arbitrary  surface  recombination 
velocity,  the  Nyquist  contribution  to  the  noise 
can  be  calculated  from  equation  (234) . 

An  accurate  approximation  to  the  excess  noise 
contribution  can  be  calculated  from  Eqs.  (243) 
and  (229)  for  the  limiting  cases  of  small  and 
large  values  of  s,  respectively.  For  large 
surface  recombination  velocity  a  triple  infin¬ 
ite  sum  is  involved  and  the  sum  appears  to 
converge  rapidly.  For  the  case  of  small  surface 
recombination  velocity  the  single  term  given  in 
Eq.  (245)  should  be  a  good  approximation.  The 
two  dimensional  solution  should  be  useful  for 
comparing  our  theoretical  results  with  experi¬ 
ment. 
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CHAPTER  IX 

SUMMARY  AND  CONCLUSIONS 


9.1  Summary  of  Procedure 

The  noise  In  a  semiconductor  having  a 
p-n  junction  results  from  the  fluctuations  of 
the  minority  carrier  density.  In  the  p-n 
junction  there  are  two  different  types  of 
minority  carriers:  holes  m  n-type  material 
and  electrons  in  p-type  material.  The  two 
types  of  minority  carriers  behave  similarly 
and  their  current  spectra  are  independent  so 
we  discuss  only  the  spectrum  of  the  hole 
density  fluctuations. 

The  solution  of  the  noise  current 
spectrum  in  a  p-n  junction  Is  complicated  if 
only  the  Kolmogorov-Fokker-Planck  (KFP)  equa¬ 
tions  are  considered.  However,  with  the  mod¬ 
ified  Langevin  approach  due  to  Petrltz2  the 
problem  is  tractable.  ThiB  method  uses  the 
KFP  equations  locally  to  derive  the  noise 
sources  and  then  uses  the  generalized  Lange¬ 
vin  equation  to  transfer  this  noise  to  the 
p-n  Junction  terminal.  In  order  to  apply 
this  method,  the  inhomogeneous  semiconductor 
equations  are  determined  using  the  transmis¬ 
sion  line-semiconductor  analogy.  From  these 
inhomogeneous  equations  the  scalar  differen¬ 
tial  equation  for  the  hole  density  and  its 
associated  scalar  Green's  function  are  found. 
The  bulk  recombination  current  spectrum  can 
now  be  solved  since  the  recombination 
process  has  only  a  scalar  source. 
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Tne  diffusion  noise  source  has  the 
characteristics  of  a  vector  source  and  the 
scalar  Green's  function  is  not  directly  applic¬ 
able.  In  solving  the  two  first  order  semicon¬ 
ductor  equations  for  the  hole  density,  the 
diffusion  source  is  changed  and  cannot  be 
determined  from  stochastic  process  theory.  We 
use  two  methods  to  remove  this  difficulty,  1, 
The  vector  current  density  equation  and  its 
associated  tensor  Green's  function  are  investi¬ 
gated.  2.  The  scalar  hole  density  equation  is 
transformed  so  that  the  diffusion  sources  can 
be  calculated  from  stochastic  process  theory. 

To  carry  out  the  first  approach,  an 
explicit  expression  for  the  tensor  Green's 
function  is  required.  A  complete  set  of  ortho¬ 
gonal  vector  eigenfunctions  is  necessary,  but 
is  not  found  for  an  arbitrary  surface  recombin¬ 
ation  velocity  on  the  transverse  surfaces.  The 
desired  set  is  found  for  the  special  case  of 
infinite  surface  recombination  velocity  on  the 
transverse  surfaces.  The  tensor  Green's  func¬ 
tion  is  determined  explicitly  and  the  diffusion 
current  spectrum  is  derived.  Adding  the  recom¬ 
bination  and  diffusion  noise  spectra,  the  total 
current  spectrum  for  infinite  surface  recombina¬ 
tion  velocity  is  determined.  This  spectrum 
consists  of  two  parts:  the  thermal  equilibrium 
Nyquist  noise  and  the  excess  noise.  Comparing 
the  N.yqulst  noise  spectrum  with  that  derived 
with  the  Nyquist  law,  we  find  that  both  spectra 
are  the  same  and  are  Infinite.  The  excess 
noise  spectrum  appears  to  be  finite  and  is  an 
upper  bound  for  the  excess  noise  spectrum  for 
semiconductors  with  large  surface  recombination 
velocity  on  the  transverse  surfaces. 

Then  to  investigate  the  second  method  of 
deriving  the  diffusion  noise  spectrum,  the 
scalar  hole  density  equation  is  transformed  by 
means  of  the  vector  integration  by  parts  and 
the  diffusion  source  takes  the  proper  form.  A 
gradient  of  the  scalar  Green's  function  appears 
which  gives  the  Green's  function  a  tensor  char¬ 
acteristic.  The  diffusion  noise  spectrum  is 
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round  for  arbitrary  surface  recombination 
velocity  on  the  transverse  surfaces.  Adding 
the  recombination  and  diffusion  noise  spectra, 
the  total  noise  spectrum  is  derived.  The 
Nyquist  part  of  this  qpectrum  is  compared  with 
the  Nyquist  noise  derived  with  the  Nyquist  law. 
Agreement  is  found  only  for  small  surface 
recombination  velocity.  The  excess  noise 
spectrum  appears  to  be  finite. 


9.2  Discussion  of  the  Noise  Current  Spectra 

The  noise  current  spectrum  depends  on  the 
area  transverse  to  the  infinite  direction;  it 
depends  on  the  characteristics  of  the  semicon¬ 
ductor  through  the  diffusion  constants,  the 
time  constants  and  the  surface  recombination 
velocity;  it  depends  on  the  geometric  shape  of 
the  semiconductor  through  the  spacial  harmonics. 
The  spectrum  is  flat  at  low  frequencies.  Only 
for  small  surface  recombination  velocity  and 
low  frequencies  is  there  a  simple  proportional¬ 
ity  between  the  excess  spectrum  and  the  dc 
longitudinal  current  at  the  p-n  junction. 

If  the  noise  spectrum  is  integrated  over 
the  whole  frequency  range,  an  infinite  noise 
current  results.  This  is  equivalent  to  the 
infinite  energy  from  a  black  body  which  was 
eliminated  by  Planck's  quantum  hypothesis. 

This  shows  that  our  equations  do  not  pertain  to 
very  high  frequencies  where  quantum  conditions 
become  Important. 

It  should  be  possible  to  compare  our 
results  with  experiment  for  the  limiting  cases 
of  large  and  small  values  of  surface  recom¬ 
bination  velocity.  For  three  dimensions  Eq. 
(217)  can  be  used  to  calculate  the  Nyquist 
noise  for  large  s,  and  Eq.  (181)  the  excess 
noise.  For  two  dimensions  the  corresponding 
equations  are  (254)  and  (229).  For  small  s 
Eqs.  (218)  and  (222)  and  Eqs.  (239)  arid  (243) 
are  the  Nyquist  and  excess  noise  contributions 
for  three  and  two  dimensions  respectively. 
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9.5  Conclusions 

In  this  thesis  we  have  contributed  to  the 
development  of  the  solution  of  a  class  of  noise 
problems.  The  method  due  to  Petritz3>4  has 
been  generalized;  using  Langevin's  determinis¬ 
tic  equations,  noise  sources  from  stochastic 
process  theory,  and  scalar  and  tensor  Green's 
functions  we  have  solved  the  p-n  junction  noise 
problem.  This  method  can  be  applied  to  other 
three  dimensional  noise  problems.  Furthermore, 
the  Green's  functions  derived  can  be  used  to 
solve  deterministic  semiconductor  problems. 

It  would  be  desirable  in  the  future  to 
extend  the  solution  obtained  with  the  tensor 
Green's  function  to  the  case  of  arbitrary  Bur- 
face  recombination  velocity.  A  set  of  orthogonal 
vector  eigenfunctions  is  required  for  this. 

Since  we  have  a  complete  set  of  independent, 
non-degenerate  eigenfunctions,  it  may  be  possible 
to  construct  a  set  of  orthogonal  vector  eigen¬ 
functions  by  means  of  Schmidt's  orthogonaliza- 
tion  process. With  a  set  of  orthogonal  vector 
eigenfunctions  it  may  be  possible  to  construct 
a  tensor  Green's  function  and  carry  out  the 
solution  of  the  problem. 

Such  a  solution,  in  addition  to  being  of 
interest  in  semiconductor  noise  theory,  would 
also  shed  some  light  on  the  disagreement  between 
our  stochastic  result  obtained  by  scalar  Green's 
function  and  the  result  obtained  by  Nyqulst's 
law. 
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APPENDIX  A 


DISCUSSION  OF  GREEN'S  FUNCTIONS 


A.l  Forward  and  Backward  Equations 

The  semiconductor  equations  as  functions  of 
time  are  affected  by  the  direction  of  time.  When 
time  increases,  these  equations  are  called  the  for¬ 
ward  equations;  Eqs.  (13)  and  (14)  with  variables 
separated  are 

V  V  W>  -  5  =  -  5“  V,t  I  ,  (Al) 


When  time  reverses,  the  semiconductor  equations  are 
called  the  backward  equations  and  are  written  as 


vv-j  («vt)  t  b  (A4 ) 
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A. 2  Equations  for  the  Green's  Functions 
and  Causality 

Each  inhomogeneous  semiconductor  equation,  (Al) 
to  (A4),  defines  an  equation  which  must  be  satisfied 
by  its  associated  Green's  function.  The  scalar  for¬ 
ward  equation  (Al)  defines  the  equation  for  the 
scalar  forward  Green's  function: 


c,  M  <*.  ,t  ,y -  ^  w  _  3^.1^  =  _  s^-osct-o  .  (as  ) 


The  scalar  backward  equation  (A3)  defines  the  equation 
for  the  scalar  backward  Green's  function: 


s/LGM\vO + % =-*(*-*}  S(t-Q.  (A6) 


DC 


The  vector  forward  equation  (A2)  defines  the  equation 
for  the  tensor  forward  Green's  function: 


^ T (r.tlr. . =  _  Jg(^)  gMo).  (A7) 


The  vector  backward  equation  (A4)  defines  the 
equation  for  the  vector  backward  Green's  function: 


«t-b.  (as) 
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Here  S  (r  -ft)  *  S  (*-  >0  S(z-7.^  (A9) 


and  S(t-tQ)  are  the  Dirac  delta  functions.1^  The 
symbol  is  the  idemfactor14  and  in  dyadic  notation  is 


,D=t,VE,t,*TA  ,  <Ai°) 

where  3,,  is  the  unit  vector  in  the  u-directlon. 

A  Green's  function  is  interpreted  as  the 
response  at  the  poirt  r  and  time  t  to  a  unit  impulse 
source  placed  at  r  and  t0.  Therefore  the  forward 
Green's  functions  Satisfy  the  causality  condition? 


lr.,to*0  ,  if  i-'t.  ,  (Ail) 

while  the  backward  or  adjoint  Green's  functions  satis¬ 
fy  the  different  causality  condition: 


IP  o  ,  if  t>t„ ,  (A12) 

¥  denotes  either  the  scalar  or  tensor  Green's 
functions. 


A. 3  The  Frequency  Domain 

When  the  semiconductor  equations  and  their  asso¬ 
ciated  Green's  functions  are  transformed  to  the 
frequency  domain,  the  forward  and  backward  equations 
become  the  same.  Using  the  technique  of  Section  3.2, 
the  scalar  equations  (Al)  and  (A3)  become  Eqs.  (34) 
and  (35),  while  the  vector  equations  become  Eqs.  (36) 
and  (37)*  Equations  (A5)  and  (A6),  which  the  scalar 
Green's  functions  satisfy,  transform  to 


vl(j"KtG  =- 6(r-ft)  ,  (A13) 

while  Eqn3.  (Ay)  and  (A8)  which  the  tensor  Green's 
functions  satisfy  transform  to 
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vS-r-tfT  — St'C  'c.'s  .  (A14) 

The  physical  concept  involved  in  the  frequency 
domain  is  that  the  sources  vary  with  a  steady-state 
sinusolsa.1  frequency.  A  steady-state  source  is  iso¬ 
tropic  in  time  and  the  forward  and  backward  equations 
are  the  same.  In  the  frequency  domain  for  a  linear 
medium  the  causality  condition  for  the  forward  Green's 
functions  is 

—  oo  4  1^(0,  ao  ,  -co  1 1  i  co  ,  (A15) 

and  likewise  for  the  backward  Green's  functions 


-col  »  -oo^-fcloo  .  (Alb) 


A. 4  Reciprocity  Relations 

Green's  functions  satisfy  a  reciprocal  relation.1^ 
The  scalar  Green's  function  for  a  single  frequency 
satisfies  the  relation 


.  (A17) 

Thl3  equation  states  that  with  a  given  harmonic 
excitation,  interchanging  the  source  and  the  observa¬ 
tion  point  does  not  affect  the  behavior  of  the  system. 

The  tensor  Oreen's  function  satisfies  a  recipro¬ 
city  relation  which  is  now  derived  (see  Morse  and 
Feshbach,  4  pp.  877-885).  Equation  (A14)  can  be 
written  as  the  operator  equation 


S (r-rft  (A18) 

where  the  differential  operator  (Lq  Is  defined  as 
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QB=  VV-  ~Kl  .  (A19) 


This  differential  operator  can  be  represented  as  the 
continuous  matrix  CL(r(ri).  With  this  matrix,  Eq. 
(A18)  becomes 


&  6H  rtV  P  C»lln,V =  -  J  & Cr-0  .  ( A20 ) 


Interchanging  rows  and  columns  and  reversing  the 
positions  of  the  factors  of  the  product  gives  the 
adjoint  of  Eq.  (A20) : 


Pfclr.)*  Cl (rl r.'j = ~J  &(r-r„)  .  (A21) 

The  tensor  Green’s  function  for  the  adjoint  operator 
31  is  the  adjoint  of  the  tensor  Green's  function 
for  CL  : 

r«=  r(r.\c)  .  (A22) 


This  is  the  reciprocity  condition  which  the  tensor 
Green's  function  must  obey. 
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APPENDIX  B 


EXCESS  CHARGE  DENSITY  AND  p-n  JUNCTION 
INPUT  ADMITTANCE 


B.l  Definition  of  Terms 

We  use  the  scalar  Green's  function,  Chapter  IV, 
to  solve  two  problems:  the  excess  hole  density 
resulting  from  a  dc  voltage  applied  to  the  plane  x»0; 
and  the  input  admittance  of  the  p-n  junction  to  an 
alternating  voltage  applied  at  x»0. 

In  the  excess  hole  density  problem  we  investi¬ 
gate  the  hole  density  at  the  x«0  plane,.  When  no 
voltage  is  applied  at  the  xsO  plane,  the  thermal 
equilibrium  hole  density  pn  exists  throughout  the 
semiconductor.  After  applying  the  voltage  V(0)  at 
the  xaQ  plane,  the  hole  density  at  xeO  increases  to 
Pt(°) • 


h  {0)  =  In  *•  P(°)  ~  +  fcieXpCfVfoj/AT)  ~'J  (Bl) 

where  q  is  the  electronic  charge,  k  is  Boltzmann's 
constant,  and  T  is  the  absolute  temperature.  The 
symbol  p  denotes  the  excess  hole  density. 

Putting  V(0)  equal  to  the  sum  of  a  large  dc 
component  VQ  and  a  very  small  ac  component  v-^, 

V(0)  =  VQ  ♦  V;L  ,  (B2) 

the  exponential  in  equation  (Bl)  is  expanded  in  a 
power  series  in  Vj_.  Keeping  the  first  two  terms,  the 
hole  density  is 

f>(o)  ft  <*) 

-H  */,  (ft/*  V,  gW'fafiffH/rrg  (w) 
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The  excess  hole  density  Is  made  up  of  two  parts: 

>  the  excess  hole  density  from  the  dc  voltage, 
VDj  and  p^(0) ,  the  excess  hole  density  from  the  ac 
voltage  v1.  Therefore, 


p  (0)  -pb( ( 0)  , 

(B4) 

fa(°)  3  fa  (flL/ATi-l], 

(B5) 

/»  (o)  =  /&r  > 

(B6) 

K>  ‘  fal^V^T)- 

(B7) 

The  average  value  of  the  total  hole  density  at 
the  x=»0  plane,  <p-t CO) > ,  is  found  by  taking  the  ensemble 
average  of  pt(0)  .  Since  the  average  value  of  the 
alternating  component  p^o)  is  zero,  we  get 


(B8) 

and  define  \  fa  ( o )fa  ~  fan 

(B9) 

and  </$  (0)}e  -  fa  (o)  , 

(BIO) 

where  the  subscripts  N  and  E  denote  the  thermal  equil¬ 
ibrium  and  the  excess  part  of  the  average  hole  density, 
respectively. 

When  a  voltage  Is  applied  at  the  x»0  plane,  an 
excess  hole  density  also  appears  throughout  the  whole 
semiconductor  and  Eqs.  (B8)  to  (BIO)  become 


(Bll) 

(B12) 

(Bl?) 
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These  equations  refer  to  finite  surface  recombination 
velocity  s  on  the  transverse  surfaces.  For  the  case 
of  infinite  s  we  put  oa  supersci’ipts  on  the  variables. 


B.2  Expass.  Hole  Density  and  Jnpu£;  Admittance 

Let  us  evaluate  the  excess  hole  density  in  the 
semiconductor  with  the  p-n  junction.  Assuming  that 
all  sources  are  zero  except  on  the  xaO  plane,  Eq.  (82) 
reduces  to 


p(r0)  (bG/3>Xo)  (Bl4) 

g  C  ifi 

P°°(r)~jo  I  />%)(*  $/* X*  >  d  V  •  (B15 ) 

where  no  superscript  and  the  &>  superscript  stand  for 
finite  and  infinite  s  cases  respectively. 

For  the  finite  s  case  (»0/ixo)  Is  calculated 
from  Eq,  (109).  When  xQ  goes  to  zero,  we  get 


*  £  Fmw  (r)  expl (Bl6) 

Since  p(rQ)  is  constant  over  the  \»0  plane,  and  Eq. 
(B5)  defines  Pb(o),  the  excess  dc  hole  density  is 


$6, 

3*. 


(c«^Xc.os^)fi^r/fc>(-V'  ^K'nrn),  (B17 ) 


where 


and  Ln,n 


A  =  4  tc,, 

is  given  by  Eq.  (ill).  This  value  of 


(B18) 

(B19) 

(B20) 

the 
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excess  hole  density  agrees  with  that  given  by 
Shockley. ' 

The  Infinite  s  case  is  similar  to  the  finite  s 
case,  and  we  get 

*  /6  rr'lfi<o)£  (mn)  (sin  irnCty  exp(-K~,^^)(b2l) 

0 


^(D?)  f  +-(rtnC')  .  (B22) 

K  nn 

The  input  admittance  to  the  p-n  junction  is 
calculated  by  finding  the  total  current  in  the  x 
direction  at  the  xsO  plane  and  dividing  by  the  ac 
voltage  V}  applied  at  this  plane.  The  current  density 
results  from  the  homogeneous  form  of  Eq.  (14)  and  from 
Eqs.  (B14)  and  (B15).  To  get  the  total  current  we 
integrate  over  the  observation  coordinates  at  the  x*0 
plane.  With  a  constant  source  function  p-,(o),  Eq, 

(B6),  at  the  x»0  plane,  the  current  in  the  x  direction 
is 


where  jjl  is  the  mobility 

(B25) 

and  the  symbols  p^,  Vj ,  G,  and  G 00  are  given  by  Eqs. 
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(B7),  (B2),  (109)  and  (114)  respectively;  the  symbols 
q,  k,  and  T  are  defined  after  Eq.  (Bl);  and  D  is 
defined  after  Eq.  (12). 

For  finite  s  the  expression  for  d*V**3<“  U.=o  is 


^57*1  KF^n^F^Cr)  KKmW  Lmv  (B26) 

M/J 

i  sc 

Putting  it  into  Eq.  (B22)  and  defining  the  input 
admittance  at  the  x«0  plane  as 


we  get 


Y| 


*  Iv, 


-« 


»-o 


X-ro  “  ^  h  % fe  ^  *«,/7 


(B27) 

(B28) 


where  A,  Jf  ,  Ki/mn ,  m.  ,  and  ph'  are  defined  by 

Eqs.  (B20),n^Bl^  (ltI?,TU2),  (B2§)  and  (B7) 
respectively. 

To  separate  Y  into  its  real  and  imaginary  partB, 


K, 


kmn 


is  so  separated.  Let 


•k 


(B29) 


where  H  is  any  complex  quantity  and  the  subscripts  r 
and  i  stand  for  real  and  imaginary  partB  respectively; 
then 


^Kmn  -44{[(Hf  *  W*  )  \  Hr  )*  i  (.(  Hr  Hr  f  j. 


.i'll 


n  . 


Using  Eq.  (B30),  Eq.  (B28)  becomes 

n..  (W&,  K‘mf 


(B30) 


(B51) 


i 
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where  Kiw-p  is  defined  by  Eq.  (B18)  and  the  other 
symbols  are  defined  after  Eq.  (B28).  This  is  the 
expression  for  admittance  which  was  derived  by 
Shockley .7 

In  the  infinite  s  case  the  input  admittance  is 


A.»  =“54  A7  ?  A »  v-  (B52) 

With,  the  real  and  imaginary  parts  shown  explicitly, 

Eq.  (BJ2)  becomes 


where ,  KjJL.  and  are  given  by  Eqs.  (B25), 

(117)  and  (B22)  respectively  and 


a  8C. 


(B34) 
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